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NONLINEAR  THOMSON  SCATTERING  OF  INTENSE  LASER  PULSES 

FROM  BEAMS  AND  PLASMAS 

I.  INTRODUCTION 

The  development  of  a  compact  source  ble,  neeir  monochromatic,  well  colh- 

mated,  short  pulse  x-rays  would  have  profouL  ide  ranging  appUcations  in  a  munber 

of  areas.  These  areas  include  x-ray  spectroscopy,  micToscopy  and  radiography,  medical 
and  biological  imaging,  x-ray  analysis  of  ultrafast  processes,  and  x-ray  holography.  One 
method  for  producing  such  an  x-ray  beam  is  by  the  nonlinear  Thomson  scattering  of  intense 
laser  pulses  from  electron  beams  and  plasmas  [1-9].  Cmrent  methods  of  x-ray  production 
include  third  generation  synchrotron  somces,  which  are  based  on  high  energy  electron 
storage  rings  and  imdulator  magnetic  fields  [10-17].  Alternatively,  x-rays  can  be  produced 
by  a  laser  synchrotron  source  (LSS),  based  on  nonlinear  Thomson  scattering,  in  which  the 
magnetic  tmdulator  is  replaced  by  ultrahigh  intensity  laser  pulses  and  the  electron  storage 
ring  is  replaced  by  a  compact  electron  accelerator  of  substantially  lower  energy  or  by  a 
stationary  plasma  [5-7].  The  compactness  of  the  LSS  makes  it  an  attractive  alternative, 
particularly  at  high  x-ray  energies  (>  10  keV),  where  conventional  synchrotrons  require 
very  high  energy  (>  5  GeV)  storage  rings.  To  generate  high  peak  fluxes  of  x-rays  in  an 
LSS,  iiltra-intense  laser  pulses  are  necessary.  Recent  advances  in  compact,  solid-state, 
short  pulse  lasers  based  on  the  method  of  chirped-pulse  amplification  [18-20],  provide  the 
technology  for  generating  the  ultrahigh  laser  intensities  required  by  an  LSS. 

In  the  following,  a  comprehensive  theory  is  developed  to  describe  the  nonlinear  Thom¬ 
son  scattering  of  intense  laser  fields  from  beams  and  plasmas.  This  theory  is  vahd  for 
hnearly  or  circularly  polarized  incident  laser  fields  of  arbitrary  intensities  and  for  electrons 
of  arbitrary  energies.  Explicit  expressions  for  the  intensity  distributions  of  the  scattered 
radiation  are  calculated  and  numerically  evaluated.  The  effects  of  the  space-charge  electro¬ 
static  potential  are  included  self-consistently  and  non-ideal  effects,  such  as  electron  energy 
spread  and  beam  emittance,  axe  discussed.  These  results  are  then  applied  to  possible  LSS 
configurations. 

An  LSS  [5-7],  using  either  an  electron  beam  or  a  plasma,  potentially  has  a  number  of 
Manuscript  approved  July  2,  1993. 
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attractive  features:  (i)  tunable  and  near-monochromatic  x-rays  can  be  obtained  over  the 
entire  x-ray  spectrum  (from  ultraviolet  to  gamma-rays),  (ii)  the  x-rays  can  be  produced 
in  ultrashort  pulses  (~  1  ps),  (iii)  a  much  lower  electron  beam  energy  (~  300  times  less) 
is  needed  to  produce  a  given  photon  energy  than  in  conventional  synchrotrons,  (iv)  the 
device  can  be  compact  and  inexpensive  compared  to  conventional  synchrotrons,  (v)  much 
higher  energy  photons  (  >  30  keV)  can  be  produced  than  in  conventional  synchrotrons,  (vi) 
the  bandwidth  can  be  small  (~  1%)  and  is  not  limited  by  the  length  of  the  undulator  as  in 
conventional  synchrotrons,  (vii)  consequently,  narrow  bandwidth  x-rays  can  be  obtained 
with  long  coherence  lengths,  (viii)  the  x-ray  polarization  is  easily  adjusted  by  changing  the 
incident  laser  polarization,  and  (ix)  high  peak  photon  fltix  and  brightness  can  be  obtained 
using  current  technology.  The  capability  of  the  LSS  in  yielding  high  average  fluxes  and 
brightnesses  is  currently  hmited  by  the  repetition  rates  of  high  intensity  laser  systems. 

An  important  parameter  in  the  discussion  of  LSS  radiation/Thomson  scattering  is  the 
dimensionless  laser  strength  parameter,  oq,  which  is  analogous  to  the  undulator  strength 
parameter,  K,  frequently  used  in  conventional  synchrotron  radiation  hterature.  The  laser 
strength  parameter  is  the  normalized  amplitude  of  the  vector  potential  of  the  incident  laser 
fleld,  ao  =  eAo/meC^,  and  is  related  to  the  laser  intensity,  Jo,  and  power,  Po,  by 

ao  =  0.85  X  10~®Ao  j^W/cm^j  (1) 

and  Po[GW]  =  21.5(aoro/Ao)^,  where  Ao  is  the  wavelength  and  ro  is  the  spot  size  of 
the  laser  (a  Gaussian  transverse  profile  is  assumed).  When  oo  1,  Thomson  scattering 
occurs  in  the  linear  regime  and  radiation  is  generated  at  the  fundamental  frequency,  u;  —  ui. 
When  oo  ^  1,  Thomson  scattering  occurs  in  the  nonlinear  regime  and  radiation  is  generated 
at  harmonics  in  addition  to  the  fundamental,  i.e.,  u>  =  ti;„  =  nivi,  where  n  =  1,2,3...  is 
the  harmonic  number.  Compact  laser  systems  based  on  chirped-pulse  amplification  can 
deliver  modest  energy  (  >  10  J),  ultrashort  (  <  1  ps)  laser  pulses  at  ultrahigh  powers  (  >  10 
TW)  and  intensities  (  >  10^®  W/cm^).  For  Ao  ~  1  /im,  oo  ^  1  requires  Iq  >  10^®  W/cm^. 
Hence,  laser  systems  which  can  be  used  to  experimentally  explore  Thomson  scattering 
in  the  nonhnear  regime  cxirrently  exist.  Furthermore,  these  powers  and  intensities  are 
sufl^cient  to  produce  ultrashort  LSS  x-ray  pulses  with  high  peak  fluxes  and  brightnesses. 
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In  the  LSS,  two  avenues  exist  for  generating  short  wavelength  radiation.  The  first  is  to 
exploit  the  relativistic  doppler  factor  which  arises  from  backscattering  laser  radiation  from 
a  counterstreaming  relativistic  electron  beam.  In  this  case,  the  wavelength  of  the  funda¬ 
mental  (n  =  1)  backscattered  radiation  along  the  axis  is  given  by  A  =  AqT^/  [(1  -I-  /5o)7o]^» 
where  70  =  {1  —  is  the  initial  relativistic  factor  of  the  electron  beam  (prior  to 

the  laser  interaction),  /3o  =  uq/c  is  the  initial  normalized  electron  velocity  and  7x  = 
(1  -I-  ao/2)^^^-  Hence,  for  70  >  1  and  Cq  •C  1,  A  ~  Ao/47o  and  extremely  short  wave¬ 
length  radiation  can  be  generated.  In  practical  units,  the  photon  energy,  Ep  =  fuD,  and 
wavelength,  A,  of  the  fundamental  backscattered  radiation  are  given  by 


0.019£||MeVl 
(l  +  ag/2)Aol,*m]’ 

A|Al=650A„w'^pg, 


(2a) 

(21.) 


where  Eb  is  the  electron  beam  energy  and  7o  1  has  been  assumed.  For  a  conven¬ 
tional  synchrotron  source  [9-16]  using  a  undulator  magnet,  A  =  Au/27o,  or  Fp[keV]  = 
0.95JS^[GeV]/At,[cm]  and  A[A]  =  13.0Au[cm]/JS5[GeV],  where  A„  is  the  imdulator  magnet 
wavelength  and  I  and  7o  1  have  been  assumed.  Since  the  laser  wavelength  in  the 
LSS  (Ao  1  /xm)  is  more  than  four  orders  of  magnitude  shorter  than  the  wavelength  of  a 
conventional  undulator  magnet  (A^,  >  4  cm),  a  much  lower  energy  electron  beam  (~  300 
times  less)  can  be  used  in  the  LSS  to  produce  a  given  photon  energy.  Hence,  compared 
to  a  conventional  storage-ring  based  synchrotron,  the  LSS  can  be  a  compact,  inexpensive 
device,  particularly  at  high  photon  energies  {Ep  >  10  keV).  As  an  example,  consider  syn¬ 
chrotron  sources  producing  30  keV  photons  (A  =  0.40  A),  assuming  Oq  -C  1  and  1. 

In  a  conventional  synchrotron  using  a  A,,  =  4  cm  undulator  period,  electron  beam  energies 
of  Eb  >  12  GeV  are  needed.  In  the  LSS  using  a  Ao  =  1  /im  laser,  Eb  =  40  MeV,  which  is 
typical  of  the  energies  available  from  compact  accelerators,  such  as  rf  linacs  or  betatrons. 

The  second  avenue  to  short  wavelengths  is  to  exploit  the  harmonic  frequency  upshift 
factor,  A  =  Ai/n,  where  Ai  is  the  wavelength  of  the  fundamental.  For  Oq  »  1,  numerous 
harmonics  are  generated.  The  result  is  a  near-continuum  of  scattered  radiation  with 
harmonics  extending  out  to  some  critical  harmonic  number,  nc  ~  Cq,  beyond  which  the 
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intensity  of  the  scattered  radiation  rapidly  decreases.  Hence,  an  ultra-intense  laser  incident 
on  a  stationary  plasma  (70  =  1)  can  generate  short  wavelength  radiation,  A  =  Xq/ti.  The 
critical  photon  energy  for  a  plasma-based  LSS  is  given  by 

Ep[eV]  =  1.24nc/Ao[/ini],  (3) 

where  ric  ~  Oq.  Assuming  laser  technology  limits  oq  ;$  10  and  Aq  ~  1  ^m  implies  that  the 
scattered  radiation  is  limited  to  A  >  10  A  and  Ep  <1  keV.  Hence,  a  plasma-based  LSS  is 
limited  by  present  laser  technology  to  the  soft  to  medium  x-ray  regime. 

Tunability  of  the  LSS  radiation  can  be  achieved  by  adjusting  either  the  electron  energy 
or  the  laser  intensity,  as  indicated  by  Eqs.  (2)  and  (3).  Neglecting  thermal  effects,  it  can  be 
shown  that  the  linewidth  of  the  scattered  radiation  for  a  particular  n  harmonic  of  frequency 
Un  is  given  by  =  1/nNo,  where  Nq  is  the  number  of  laser  periods  with  which  the 

electron  interacts.  In  principle,  since  Nq  is  typically  large  {Nq  >  300),  narrow  linewidth 
x-rays  can  be  generated.  In  practice,  the  linewidth  will  be  limited  by  thermal  effects. 
For  example,  the  normalized  energy  spread  associated  with  an  electron  beam,  AE/Ef,, 
limits  the  linewidth  to  2AE/Eb.  An  additional  advantage  of  generating  LSS 

radiation  using  an  electron  beam  is  that  the  scattered  radiation  is  weU  collimated  about  the 
backscattered  direction  (i.e.,  the  direction  of  the  electron  beam).  For  an  electron  beam  with 
7o  ^  1  and  cq  <  1,  the  backscattered  radiation  with  linewidth  Au;/u;  ~  I/Nq  is  confined 
to  a  radiation  cone  of  half-angle  6  ~  l/(7o\/^).  For  a  plasma  with  oq  >  1,  the  radiation 
is  scattered  over  a  much  larger  angle.  When  ao  1,  numerous  harmonics  are  generated, 
and  tunability  is  achieved  by  filtering  the  scattered  radiation.  An  additional  advantage 
in  using  a  plasma  is  that  very  high  electron  densities  can  be  achieved  in  comparison  to 
densities  obtainable  in  electron  beams.  The  scattered  power,  as  well  as  photon  flux  and 
brightness,  scale  linearly  with  density,  hence,  the  use  of  high  electron  densities  is  favored. 

Thomson  scattering  theory  is  a  classical  description  which  is  valid  provided  the  scat¬ 
tered  photon  energy  is  small  compared  to  the  electron  energy,  i.e.,  fiu;  <  7oTneC^.  For  a 
plasma,  this  implies  photon  energies  less  than  500  keV.  For  an  electron  beam  with  70  1, 

Ao  =  1  /im  and  oo  <  1,  this  implies  70  <  10®,  i.e.,  electron  beam  energies  less  than  50 
GeV.  Nonlinear  Thomson  scattering  of  intense  radiation  from  a  single  electron  initially 
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at  rest  was  examined  analytically  in  considerable  detail  in  the  classic  work  of  Sarachik 
and  Schappert  [1].  (This  work  was  recently  reexamined  by  Castillo-Herrera  and  Johnston 
[9].)  However,  the  important  effects  of  the  space  charge  potential  [2,21],  which  arises  in 
high  density  plasmas,  was  neglected  and  scattering  from  electron  beams  was  not  discussed. 
Waltz  and  Manley  [2]  also  discussed  Thomson  scattering  from  plasmas  and  pointed  out 
that  the  space  charge  potential  was  important  in  preventing  the  drift  of  electrons  in  the 
direction  of  the  incident  laser.  However,  explicit  expressions  for  the  scattered  intensity  dis¬ 
tribution  for  arbitrary  ao  were  not  calculated  and  scattering  from  electron  beeims  was  not 
considered.  Many  authors  [10-17]  have  analyzed  the  production  of  synchrotron  radiation 
in  the  interaction  of  relativistic  electron  beams  vdth  static  magnetic  imdulator  and  wiggler 
fields,  a  process  which  is  somewhat  similar  to  Thomson  scattering.  These  analyses  require 
that  if/7o  1  (analogous  to  ao/70  ^  1)>  2^  assumption  which  need  not  be  made  in  the 
analysis  of  nonlinear  Thomson  scattering.  In  this  paper,  nonlinear  Thomson  scattering 
of  intense  laser  fields  from  electron  beams  and  from  plasmas  is  examined  analytically  and 
numerically.  This  analysis  is  valid  for  linearly  and  circularly  polarized  incident  laser  fields 
of  arbitrary  intensities  and  for  electron  beams  of  arbitrary  energies  (up  to  the  limits  of 
classical  theory).  The  effects  of  the  space-charge  potential  are  included  self-consistently 
and  various  non-ideal  effects,  such  as  electron  energy  spread,  are  discussed. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Sect.  H,  the  orbits  of  electrons 
in  intense  laser  fields,  both  linearly  and  circularly  polarized,  are  calculated  including  the 
effects  of  the  self-consistent  electrostatic  potential.  Explicit  expressions  for  the  scattered 
intensity  distributions  are  derived  in  Sect.  III.  These  are  genered  expressions,  valid  for 
electron  beams  and  plasmas,  and  for  arbitrary  laser  intensities.  Properties  of  the  scattered 
radiation  are  examined  in  Sect.  IV,  including  a  calculation  of  the  total  power  radiated  from 
an  electron  beam  or  a  plasma,  an  examination  of  the  resonance  function  and  the  behavior  of 
the  radiation  spectra  in  the  ultra-intense  regime,  i.e.,  Oq  »  1.  Various  non-ideal  effects  are 
discussed  in  Sect.  V,  including  the  effects  of  electron  energy  spread,  electron  beam  energy 
loss,  ponderomotive  density  depletion  and  plasma  dispersion.  These  results  are  apphed  to 
possible  LSS  configurations  in  Sect.  VI,  and  specific  examples  of  an  electron-beam  LSS 
and  a  plasma  LSS  are  presented.  Section  VII  is  the  conclusion. 
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II.  ELECTRON  MOTION  IN  INTENSE  LASER  FIELDS 

The  laser  field  and  space  charge  field  of  the  electrons  can  be  represented  using  the 
normahzed  vector  and  scalar  potentials,  a  =  eA/meC^  and  $  =  respectively, 

where  me  is  the  electron  mass  and  c  is  the  magnitude  of  the  electron  charge.  In  the 
Coulomb  gauge,  V  •  a  =  0  implies  =  0  in  one-dimension  (ID).  Then,  oj.  represents  the 
laser  field  and  $  represents  the  space-charge  field  of  the  plcisma.  The  normalized  vector 
potential  of  a  laser  of  arbitrary  polarization  is  represented  by 


a  =  (ao/\/2)  ^(1  +  cos  koT}  Cx  +  (1  -  sin  koV  Byj  ,  (4) 

where  fco  =  Stt/Aq  is  the  wavenumber  of  the  laser  field,  tj  =  z  +  ct,  6p  =  1  for  Unear 
polarization  and  6p  =  0  for  circular  polarization.  Using  this  representation,  (a^)s  = 
ao/2  for  both  linear  and  circular  polarizations,  where  the  subscript  s  signifies  the  slow 
component  (an  averaging  over  the  laser  wavelength).  Hence,  the  average  laser  power 
Pq  ~  (a^)s  is  constant  for  a  given  value  of  uq,  independent  of  polzurization,  i.e.,  Po[GW]  = 
21.5(aoro/Ao)^,  assuming  a  Gaussian  transverse  profile  of  the  form  ja]  ~  exp(— r^/rp).  In 
the  following,  the  laser  field  is  assumed  to  be  moving  to  the  left  {-z  direction)  and  the 
electrons  are  initially  (prior  to  the  interaction  with  the  laser  field)  moving  to  the  right  (+z 
direction)  with  an  initial  axial  velocity  u*  =  vq  (see  Fig.  1). 

The  electron  motion  in  the  fields  a  and  $  is  governed  by  the  relativistic  Lorentz 
equation,  which  may  be  written  in  the  form 

1  d  19  _ 

u  =  V$+ a-/3x  (Vxa),  (5) 

c  at  c  at 

where  /3  =  v/c  is  the  normalized  electron  velocity,  u  =  p/meC  =  7/3  is  the  normalized 
electron  momenttun,  and  7  =  (1+w^)^/^  =  is  the  relativistic  factor.  Assuming 

that  the  laser  field,  ax,  and  hence  the  quantities  4,  /3,  u,  and  7,  are  functions  only  of  the 
variable  77  =  z  +  ct,  Eq.  (5)  impUes  the  existence  of  two  constants  of  the  motion  [21,22], 

^(ux-ax)  =  0,  (6a) 

4-{l  +  -  ^)  =  0.  {6b) 
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Equation  (6a)  is  conservation  of  canonical  transverse  momentum  in  ID,  and  Eq.  (6b)  can 
be  interpreted  as  conservation  of  energy  in  the  wave  frame.  Equations  (6a)  and  (6b)  can 
be  integrated  to  give  [21,22] 


Ux  =  ax, 

(7a) 

7  +  u,  -  $  =  7o(l  +  0q), 

(7b) 

where,  prior  to  the  laser  interaction  (ax  =  0),  ux  =  ^  =  0,  7  =  70  and  Uz  =  7o/?o  have 
been  assumed.  The  two  constants  of  the  motion,  Eqs.  (7a)  and  (7b),  completely  describe 

the  nonlinear  motion  of  electrons  in  the  potentials  a  and  They  allow  the  electron  motion 

to  be  specified  solely  in  terms  of  the  fields,  i.e.. 

^  hg  -  (1  +  a^) 
^^-hi+(i+a^r 

(8a) 

7  =  (^0  ■i"  1  d"  a^)/2/io, 

(86) 

Px  =  a±/7, 

(8c) 

where  ho  =  70(1  +  0o)  +  $. 

The  self-consistent  space-charge  potential  of  the  electrons, 

the  continuity  equation  and  Poisson’s  equation, 

can  be  determined  using 

-^ne  +  V(ne/3)  =  0, 
c  eft 

(9a) 

V2$  =  fcJ(ne/no-l), 

(96) 

where  Ue  is  the  electron  density,  kp  =  Wp/c,  Wp  =  {4Tre^no/mey^^  is  the  plasma  frequency 
and  no  is  the  ambient  density.  Equation  (9b)  assumes  that  the  initial  equilibrium  (prior 
to  the  laser  pxilse)  space-charge  potential,  is  negligible.  For  a  plasma,  a  neutralizing 
background  of  stationary  ions  is  assumed,  i.e.,  =  0.  For  a  long,  uniform  electron 

beam  of  radius  rj,,  <  fcJr^/4  =  Ub,  where  Ub  =  h/Ibo  is  the  Budker  parameter,  h  is 

the  beam  current,  and  IwlkA]  =  17 0^.  Since  i/j,  1  for  beams  of  interest,  can  be 
neglected.  Assuming  rie  =  rieiil),  Eq.  (9a)  implies  [21,22] 

^Mi  +  a)]  =  o,  (10) 
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hence,  n*  =  no(l  +  /3o)/U  +  0z)-  Substituting  this  result  into  Eq.  (9b)  and  using  Eq.  (8a) 
give  [21,22] 


2  [(1  +  'J')2 


(11) 


where  ^  -  4/7o(l  +  6o)  and  kp  =  A:p/7o^^(l  +  ^o)- 

Equation  (11)  describes  the  self-consistent  electrostatic  potential  induced  by  the  in¬ 
teraction  of  the  laser  field.  The  solution  for  is,  in  general,  highly  nonlinear.  Simple 
solutions  can  be  obtained  in  two  limits  in  which  the  characteristic  temporal  variation  of 
the  laser  envelope,  tl  (typically  the  laser  rise  time),  is  compared  to  an  effective  plasma 
period,  {ckp)~^.  In  the  short-pulse  limit,  r^,  C  (cfcp)“^,  Eq.  (11)  imphes  |'i'|  ■C  1  provided 
og  <  2/cr£,fcp,  where  uq  is  the  amphtude  of  the  laser  pulse,  e.g.,  a  =  oq  cos  In  the  long 
pulse  limit,  tl  »  {ckp)~^,  the  left  side  of  Eq.  (11)  can  be  neglected  and  it  can  be  shown 
that  ^  ~  (1  +  —  1,  where  the  subscript  s  signifies  the  slow  part.  Throughout  the 

following,  the  quantity  (H-a^)^^  ~  (l-|-ao/2)^/2  be  approximated  as  nearly  constant, 
i.e.,  \d{a^)sldTi\  ko(a^)a,  which  implies  that  Lq  ^  Aq,  where  Lq  =  cti  is  the  length  of 
the  laser  envelope. 

For  applications  which  utilize  intense  lasers  with  pulse  lengths  rx,  ~  1  ps,  the  short- 
pulse  limit  is  relevant  to  interactions  with  electron  beams  as  long  as  the  beam  density  is 
sufficiently  low,  no/7o  ^  10^®  cm“^.  On  the  other  hand,  the  long-puke  limit  is  relevant 

to  interactions  with  stationary  (70  =  1)  plasmas  as  long  as  the  density  k  sufficiently  high, 

1 

no  2>  10^®  cm“®.  Under  these  conditions,  the  parameter  ho  =  70(1  +  A))(l  +  ’5')  k  given 
by 

[  7o(l  +  /3o),  e-beam  (short  puke), 
ho  =  s  (12) 

[  (1  -I-  Oo/2)^/2^  plasma  (long  puke). 

Notice  that  in  the  limit  of  a  low-density  plasma  with  no  -C  10^®,  |4l  1  and  Hq  ~  1. 

This  corresponds  to  the  single  particle  limit  considered  in  Ref.  [1]. 

The  electron  orbits,  T{ri)  =  x  e,  -I-  y  e,,  -I-  z  e,,  can  be  calculated  as  a  function  of  t} 
using  Eqs.  (8a-c)  and  the  relation 


1  dr 
c  dt 


=  3=(i+A)|, 


(13) 
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which  gives  dr /dr)  =  u/Ziq.  For  a  hnearly  polarized  laser  of  the  form  given  by  Eq.  (4)  with 
tfp  =  1,  the  electron  orbits  are  given  by 


Ux  =  Oo  COS  koT), 

(14a) 

Uy  —  0) 

(146) 

Uz  =  [ho  -  (1  +  OflCos^  koT))]  /2ho. 

(14c) 

Hence, 

x{t))  =  XQ  +  ri  sin  kov. 

(15a) 

II 

(15b) 

z{t))  =  zq  +  PiT)  +  zism  2koT), 

(15c) 

where  additional  terms  of  order  Xq/Lq  have  been  neglected  and 

ri  =  ao/hoko, 

(16a) 

=  -al/Shlko, 

(16b) 

/3i  =  (l-l/Mo)/2, 

(16c) 

with  Mo  =  h^/(l  +  ag/2),  i.e., 

(  7^(1  +  0of/il  +  al/2),  e-beam, 

Mo=  ^ 

(17) 

y  1,  plasma. 

Similarly,  for  a  circular  polarized  Itiser  {6p  =  0),  the  electron  orbits 

are  given  by 

itx  =  {ao/\/2)coskoV, 

(18a) 

Uy  =  (ao/\/2)sinfco^, 

(186) 

tiz=  W-(l+ai/2)]  /2ho. 

(18c) 

Hence, 

x(t/)  =  xo  +  (ri/V^)sinfco»7, 

(19a) 

yiv)  =  yo~  {r\lV2)  cos  kor), 

(196) 

z{r))  =  20  + 

(19c) 
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where,  agadn,  additional  terms  of  order  XofLo  have  been  neglected.  In  the  above  equations, 
(^JQ)  2/0.  -^o)  are  related  to  the  initial  position  of  the  electron. 

The  axial  drift  velocity  of  the  electrons,  fit ,  can  be  written  in  terms  of  the  pmameter 
fii-  Since  77  =  z  +  ct,  Eq.  (19c)  imphes  z  =  (zq  +  fiict)/{l  —  fii).  Hence, 

fit  =  fii/{l  -  fix)  =  (Mo  -  l)/(Mo  +  1)  (20) 

is  the  average  normalized  velocity  of  the  electrons  in  the  axial  direction.  Notice  that  in 
the  dense  plasma  (long  pulse)  hmit,  Mo  =  1  and  fit  =  0.  For  a  low  density  plasma  in  the 
single  particle  limit  Mo  =  (H-ao/2)~^  and  fit  —  — (ao/2)/(2 -f  ao/2).  Hence,  in  the  single 
particle  hmit,  a  single  electron  initially  at  rest  receives  a  finite  average  drift  velocity  due 
to  the  ponderomotive  force  associated  with  the  rise  of  the  incident  laser  pulse,  as  pointed 
out  in  Ref.  [l].  For  an  electron  in  a  dense  plasma  (long  pulse  limit),  fit  =  ^  and  there  is 
no  average  axial  motion  of  the  electrons  [2,21,22].  Physically,  =  0  is  achieved  through 
a  balance  between  the  ponderomotive  force  and  the  space-charge  force  set  up  during  the 
rise  of  the  laser  pulse. 

The  self-consistent  electron  density  in  the  presence  of  the  laser  field  can  be  calculated 
using  the  constant  of  motion  ne(l  +  fit)  =  no(l  +  /^o)-  This  can  be  written  in  terms  of  the 
parameter  ho  ^ 

=  wo(l  + /3o)(ho  +  1  +  a^)/2ho-  (21) 

Of  particular  interest  is  the  slow  part  (t?  averaged)  of  the  density,  Uej.  For  a  tenuous 
electron  beam  (short  pulse  hmit),  ho  =  7o(l  +  /3o)  aJid  ne*  ^  1,  assuming  hp  >  (H-ao/2). 
For  a  dense  plasma  (long  pulse  hmit),  ho  =  (1  +  ao/2)^/^  and  Ties  =  1-  However,  this  is 
not  the  case  for  a  plasma  in  the  single  particle  regime.  For  a  tenuous  plasma  in  the  short 
pulse  hmit,  ho  =  1  and  Ties  =  «o(l  +ao/4)-  regime,  the  plasma  density  is  enhanced 

due  the  ponderomotive  force  associated  with  the  rise  of  the  laser  pulse  and  the  resulting 
finite  axial  drift  motion  of  the  electrons,  fit- 

The  above  results  have  assumed  the  ID  hmit,  which  is  vahd  when  tq  »  Ao  and  when 
the  quiver  motion  is  much  greater  than  the  ponderomotive  motion.  In  three-dimensions 
(3D),  the  ponderomotive  motion,  =  u  —  a,  is  given  [23]  by  dbixldr}  =  V((/»  -  7).  The 
quasi-static  approximation  imphes  that  the  quantity  7  +  u,  —  0  —  a,  is  a  constant  of 
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the  motion,  which  is  the  3D  generahzation  of  Eq.  (7b).  For  a  plasma,  it  follows  that 
l<5“l/!a|  Apflo/ro,  whereas  for  a  relativistic  electron  beam,  |6u|/|a|  <,  I,oao/7o^o-  The 
ponderomotive  motion  can  be  neglected  when  |6ii|/|a|  1,  which  is  true  in  the  cases 

discussed  below. 
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III.  SCATTERED  RADIATION 


The  energy  spectrum  of  the  radiation  emitted  by  a  single  electron  in  an  aurbitrary  orbit 
r(t)  and  P{t)  can  be  calculated  from  the  Lienard-Wiechert  potentials  [24], 


JlL 

dudfl 


I  r 

47r2c  I  J_ 


r/2  ,2 

dt  [n  X  (n  X  ^)]  exp  —  n  •  r/c)] 
r/2 


(22) 


where  d^I/dujdfl  is  the  energy  radiated  per  frequency,  iv,  per  sohd  angle,  ft,  during  the  in¬ 
teraction  time,  T,  and  n  is  a  unit  vector  pointing  in  the  direction  of  observation.  Introduc¬ 
ing  the  spherical  coordinates  {r,6,<f))  and  unit  vectors  (er,ee,e^),  where  x  =  rsui6cos(f>, 
2/  =  r  sin  0  sin  0,  z  =  r  cos  6,  and 


Cr  =  sin  d  cos  0  Ca;  -f  sin  ^  sin  Cy  -I-  cos  6  ,  (23a) 

Cfl  =  cos  6  cos  <t>  ex  +  cos  0  sin  Cy  —  sin  0  e* ,  (236) 

Gti,  =  —sni<f)  Bx  +  cos  (f>  Cy,  (23c) 

and  by  identifying  =  n,  give 

n  X  (n  X  /3)  =  -{(3x  cos  0  cos  (/» -f  Py  cos  0  sin  sin  6)  eg 

-t-  (fix  fiy  cos  <^)  e^,  (24a) 

n  •  r  =  ar  sin  ^  cos  (^-|-ysin^sin0-l-z  cos  9.  (246) 


The  scattered  radiation  will  be  polarized  in  the  direction  of  n  x  (n  x  fi).  Hence,  I  =  -f  7,^, 
where  Jg  aind  7,^  are  the  energies  radiated  with  polarizations  in  the  eg  and  directions. 


respectively.  In  terms  of  the  independent  variable  r}  =  z  +  ct, 


d'^Ig 

dudfl 


I  f  dx  dy  dz  \  .|2 

J  dr/ cos5cos0+ ^cos^sin^— —  sin^J  exp(iV’)|  ,  (25a) 


_ 
4ir^(^  I 


dudfl 


e^u^ 


47r^c^ 


(256) 


where 


•0  =  fc  [r/  —  2(1  -f  cos  6)  —  X  sin  9  cos  0  —  y  sin  0  sin  0] , 


(26) 


k  =  u/c,  r/o  =  Lo/2,  Lq  is  the  laser  pulse  length  and  »  Ao  =  27r/fco  has  been  assumed. 
In  deriving  the  above  expressions,  the  relation  cfidt  =  {dT/dT])dT]  was  used,  where  r  =  r(T/) 
is  given  by  Eqs.  (15)  and  (19). 
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A.  Linear  Polarization 


The  electron  orbit  for  a  lineaurly  polarized  incident  laser  field  of  the  form  given  by  Eq. 
(4)  with  6p  =  1  is  given  by  Eq.  (15).  The  phase  ^  can  be  written  as 

V’  =  ■00  +  [1  -  A(1  +  cos  6)]  kq 

—  {kri  sm$cos<f>)smkoq  —  [kzi{l  +  cos ^)]  sin  2/:o^,  (27a) 

00  =  — [20(1  +  COS  d)  +  xo  sin  0 cos  0  +  sin  ^ sin  0] .  (275) 

Using  the  Bessel  identity 

00 

exp(i6sincr)  =  ^2  •7n(5)  exp(m<T),  (28) 

n=*-oo 

where  J„  are  Bessel  functions,  allows  the  phase  factor  exp  [i(0  +  £koq)]  to  be  written  as 

00 

exp  [i(0  +  £koq)]  =  ^  •7m(Q:z)*7n-2m+/(a*)  exp  [i(0o  +  fc;/)]  ,  (29) 

m,n=— 00 

where 


=  A:  [1  -  /?i(l  +  cos  9)]  —  nko,  (30a) 

Oz  =  kzi(l  +  cos  9),  (305) 

dx  =  kri  sin  9  cos  0.  (30c) 

In  order  to  evaluate  Eqs.  (25a)  and  (25b),  it  is  necessary  to  evaluate  the  integrals 

=  J  exp(i0).  (31) 

Using  the  orbits,  Eq.  (15),  along  with  the  identities  in  Eqs.  (28)  and  (29), 

L  =fcoT’ie‘^°  f;  Jmiaz)  (Jn-2m-l(dx)  +  Jn-2m+l (d*)] ,  (32a) 

m,n=— oo  '  ' 

/.=2e-Vo  ±  (^)j™(a.) 

m,n=— oo  '  ' 

•  {PlJn-2m{Ox)  +  fcoZl  [•7n-2m-2(d*)  +  Jn-2Tn+2(d*)]}  ,(326) 


and  iy  =  0,  where 


iPle  ^ _ 

dtudn 

<Pl^  _ 


ix  cos  6  cos  <l>  —  iz  sin  6 
2 

JxSin^ 


(зза) 

(ззб) 


dujdO. 

The  frequency  width  of  the  radiation  spectrum  for  a  given  harmonic  is  determined  by 
the  resonance  function  R{k,nko),  where 

(I?)'- 

This  function  is  sharply  peaked  about  the  resonant  frequency,  u;„,  given  by  fc  =  0, 


nu>o _ 

1  —  0i(l  +  COS0)' 


(35) 


The  width  of  the  spectrum,  Au,  about  cjn  is  given  by  Au/un  =  1/nNo,  where  iVo  =  Lq/Xq 
is  the  number  of  periods  of  the  laser  field  with  which  the  electron  interacts. 

Since  the  firequency  spectra  for  two  difierent  harmonics,  n  and  n\  are  suflSciently  well 
separated,  the  summations  in  Eqs.  (32a)  and  (32b)  may  be  simplified  to  yield 

/sinfct;o'\^ 

dujdSl  ^A'k'^c  \  ifc  / 

■  [Ca(l  —  sin^  6  cos^  4>)  +  sin^  9  —  CxC^  sin  29  cos  <t>\  ,  (36) 


where 

OO 

Cx  —  ^  ]  (“1)  ^0^1«An(^z)  («Ai— 2m— l(o£x)  “I"  •Ai— 2m+l (o^x)]  i  (37a) 

m=— OO 

OO 

Y.  (-ir2Jm(aJ{i8iJn-2m(ax) 

m=— OO 

+  k^Zi  \Jn—2m—2iS^x)  "I"  *Ai— 2m+2(Q^i)]  }>  (376) 


and 


_  nao(l  +  cos6) 

8hg  [1  -  /3i(H-  cos  0)]  ’ 
nao  sin  9  cos  0 
*  ^  ho  [1- A(l  +  cos^)]‘ 


(38a) 

(386) 
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In  the  deriving  the  above  expressions,  the  approximation  w  ~  was  made  in  the  argu¬ 
ments  of  the  Bessel  functions,  Qx  and  a^. 

Plots  of  the  normcilized  amplitude  of  the  scattered  intensity,  d^I  jdud^,  versus  normal¬ 
ized  frequency,  a;/47ou;o,  and  normahzed  observation  angle,  7o0,  are  shown  in  Figs.  2(a,b) 
for  the  case  of  a  linearly  polarized  laser  (iVo  =  7)  interacting  with  a  counterpropagating 
relativistic  electron  (70  =  5).  The  intensity  is  shown  in  the  plane  of  electron  motion,  </>  =  0, 
i.e.,  9  is  the  “horizontal”  observation  angle  (0  =  0  is  along  the  z-cixis,  the  axis  of  propaga¬ 
tion).  Figure  2(a)  shows  the  intensity  in  the  first  two  harmonics  for  oq  =  0.5.  Significant 
radiation  occurs  only  at  the  fundamental  (n  =  1).  The  intensity  of  the  fundamental  peaks 
on  axis  with  a  frequency  shifted  slightly  from  the  low-intensity,  Thomson  backscattered 
value  of  47oaJo,  and  is  confined  to  an  angle  6  <  I/70.  Figure  2(b)  shows  the  intensity  in 
the  first  three  harmonics  for  Oq  =  1.0.  Significant  radiation  now  occurs  in  the  harmonics  as 
well  as  the  fundamental.  Only  the  odd  haurmonics  are  finite  adong  the  axis  (6  =  0)  and  the 
frequency  shift  due  to  finite  oq  is  more  apparent.  The  angular  distribution  of  the  higher 
harmonics  is  more  extensive  than  the  fimdamental.  The  harmonic  exhibits  (n-hl)/2, 
for  n  odd,  or  n/2,  for  n  even,  intensity  maxima  as  a  function  of  6.  For  larger  values  of  oo, 
the  harmonics  dominate  the  spectrum. 

Plots  of  the  normalized  amplitude  of  the  scattered  intensity,  d'^I/dudfl,  versus  obser¬ 
vation  angle,  6,  are  shown  in  Figs.  3(a-c)  for  the  case  of  a  linearly  polarized  laser  {Nq  =  7) 
interacting  with  a  dense  plasma  electron.  The  intensity  is  shown  in  the  plane  of  electron 
motion,  <f)  =  0,  i.e.,  9  is  the  horizontal  observation  angle.  Figure  3(a)  shows  the  intensity 
in  first  three  harmonics  for  oq  =  0.5,  Fig.  3(b)  shows  the  intensity  in  first  six  harmonics 
for  Oo  =  1.0,  and  Fig.  3(c)  shows  the  intensity  in  first  twelve  h2U’monics  for  oq  =  2.0.  For 
a  dense  plasma,  there  is  no  average  axial  drift  of  the  electrons,  hence,  heu'monic  radiation 
is  scattered  over  large  angles  and  the  frequency  is  not  shifted,  i.e.,  =  tkjq-  (For  conve¬ 

nience,  the  intensity  is  plotted  only  at  the  resonant  frequencies,  u>  =  u;„.)  Only  the  odd 
harmonics  are  finite  along  the  axis  (9  =  0)  and  the  intensity  is  maximum  off-axis  for  all 
harmonics  with  n  >  1.  The  n*'*  harmonic  exhibits  (n  -I- 1)/2,  for  n  odd,  or  n/2,  for  n  even, 
intensity  maxima  as  a  function  of  9  within  the  region  0  <  9  <  7r/2. 

Backscattered  Radiation.  Of  particular  interest  is  the  radiation  backscattered  along 
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the  aixis.  In  the  backscattered  direction,  6  =  0,  only  the  odd  harmonics  axe  finite, 
the  even  harmonics  vanish.  Setting  6  =  0  in  the  above  expressions  gives,  for  the 
harmonic. 


d^In 


e=o 


c2fcoiVoMo"F„(ao)G„(u;), 


i.e., 

odd 


(39) 


where 


Fn(ao)  =  nan  [j(„-l)/2(an)  -  •/(n+l)/2(an)]” 


(40) 


is  the  harmonic  amplitude  function,  an  =  nao/4(l  +  0^/2), 


\ 

Gn{(*^) 


R{k,nko) 

Aui 


1  sin(u;  —  nMcoJo)? 
Aa;  [  (uJ  —  nMou}o)T 


(41) 


is  the  frequency  spectrum  function  and  T  =  Lo/2cMo.  The  function  G„(£i;)  is  a  resonance 
function  sharply  peaked  about  the  resonant  frequency,  a;„  =  nMouio,  with  a  width  given 
by  Auijuin  =  l/wiVoi  where  the  frequency  multiplication  factor  Mq  is  given  by  Eq.  (17). 
Furthermore,  6(a;  —  a;„)  as  iVo  — >  oo. 

The  energy  radiated  in  the  backscattered  harmonic  depends  on  the  fimction 
Fn(ao))  Eq-  (40).  For  high  harmonics,  n  >  1,  F„  becomes  significant  when  Oq  »  1. 
For  modest  power  lasers  for  which  Oq  -C  1,  only  the  fundamental,  n  =  1,  is  significant.  A 
plot  of  the  function  Fn  versus  the  parameter  (aQ/4)/(l  +  Oo/2)  is  shown  in  Fig.  4. 

B.  Circular  Polarization 


To  calculate  the  scattered  radiation  from  a  circularly  polarized  incident  laser  field 
(^p  =  0),  the  orbits  given  by  Eqs.  (18)-(19)  are  used  in  Eqs.  (25a)-(25b).  The  intensity 
distribution  can  be  written  as 


d^h 


0 


dudSl 

d^I^, 

du)dCl 


47r2c^  I  y_ 

»2,  ,2 


>70 


dr) 


VO 

VO 


kon 


cos( 


=  g  I  f 

47r2c3  I 


dr] 


^  ^  ^^^OcosikoT]  01  sin^j  exp(i'0)|  , 

exp(tV')|  . 


-  <l>) 


(42a) 

(426) 


The  phase,  V*,  is  given  by 


^  -  iPq  ■+ [1  -  /3i{l  +  COS0)]  kr)  -  (fcri/V^)sin0sin(A:oT?  -  <f>),  (43) 
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where  ■0o  is  given  by  Eq.  (27b).  Using  the  Bessel  identity,  Eq.  (28),  gives 


exp{i[V’  +  ^(A;or?- <;6)]}  =  exp  [rXt/'o  +  ^^  +  «<^)] ‘^n+<(a), 


where  d  =  {kr^/y/i)sm.6  and  where  k  is  given  by  Eq.  (30a).  This  allows  the  calciilation 
of  the  integrals  in  Eqs.  (42a)  and  (42b).  In  particular. 


/Vo 

drjexpiiip) 

■Vo 

=  ^  exp  +  n0)]  2J„(q:), 

n=— oo  '  ' 

/Vo 

drj  cos(A:o»7  —  V*)  exp{irp) 

■Vo 

=  ^  exp  [i(V>o  +  n<f))]  ^  ^  Jn(d), 

n=— oo  '  ' 

/VO 

dr)  sm{kor}  -  rp)  exp(iV') 

■VO 

=  5^  exp  [t(V»o  +  n((>)]  2iJ;(q:). 


(45a) 


(456) 


(45c) 


As  indicated  by  Eq.  (34),  the  above  expressions  imply  a  jfrequency  spectrum  centered 
about  (jj  =  (jJn,  where  w„  is  given  by  Eq.  (35),  of  width  Aw/wn  =  I/tiNq.  Since  the 
frequency  spectra  of  two  different  harmonics,  n  and  n',  are  well  sepeirated,  the  summations 
in  Eqs.  (42a)-(42b)  can  be  simplified.  Using  Eqs.  (42)  and  (45),  the  radiation  spectrum 


can  be  written  as 

d?I 

dudO, 


_  ■^e^k^  ( sinkT^oV 

”  “i  V  k  ) 


cos 6  —  cos oy^ 

'  2  £k 

Sin  d 


Mri  ^2,  ,1 

k®)  "b  2  I  ’ 


where  kovi  =  ao/ho  and  the  approximation  w  ~  u;„  has  been  made  in  the  arguments  of 
the  Bessel  functions,  i.e., 

^  n(ao/t/2)sing 
ho  [1  -  Pi{l  +  COS0)]' 
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In  the  above  expression,  the  terms  proportional  to  J„(a)  are  the  contributions  from  1$, 
and  the  terms  proportional  to  J^(a)  are  the  contributions  from 
Using  the  identities  [1] 


Y^n^Jn(nz) 

n=:l 


i2(4  +  i2) 

16{l-i2)7/2’ 


5^  Jn^(ni)  = 

n=l 


(4  +  3z2) 

16(1-22)5/2’ 


(48) 


the  summation  in  Eq.  (46)  can  be  carried  out  and  an  expression  for  dl/dft  can  be  found. 
After  integrating  over  frequency,  one  finds 


dO, 


{e'^/c)NoLJoal/hl 


32(1  -  £2)7/2  [1  -  /3i(l  +  cos^)]^ 


f  [cos0  —  /3i(l  +  cos0)]^ 
\  [1  - /3i(l +  cos0)f 


(4  +  z^)  +  (4  +  3i^) 


(49) 


where  i  =  a/n. 

Plots  of  the  normalized  amplitude  of  the  scattered  intensity,  d^I/du)dn,  versus  nor¬ 
malized  frequency,  u/4'yQUo,  and  normalized  observation  angle,  'yod,  are  shown  in  Figs.  5 
and  6  for  the  case  of  a  circularly  polarized  laser  (ao  =  1,  Ao  =  7).  Because  of  the  symmetry 
of  the  electron  orbit,  the  intensity  distribution  is  independent  of  <f>.  Figure  5  shows  the 
scattered  intensity  from  a  counterpropagating  relativistic  electron  (70  =  5)  for  the  first 
three  harmonics.  Only  the  fundamental  (n  =  1)  is  nonzero  on  axis,  where  its  intensity 
is  maximum,  and  its  frequency  is  shifted  from  the  low-intensity,  Thomson  backscattered 
value  of  47ou;o.  The  intensity  of  the  higher  harmonics  peak  off-axis  and  is  confined  to  an¬ 
gles  9  <  as  discussed  in  Sect.  IV  C  below.  Figure  6  shows  the  scattered  intensity 

from  a  electron  in  a  dense  plasma  for  the  first  six  harmonics.  For  a  dense  plasma,  there 
is  no  average  axial  drift  of  the  electrons  and  the  frequency  is  not  shifted,  i.e.,  u;„  =  rujo- 
Only  the  fundamental  is  nonzero  on  axis,  where  its  intensity  is  maximum.  For  higher 
harmonics,  the  intensity  is  maximum  in  the  transverse  direction,  =  7r/2.  As  the  intensity 
of  the  laser  pulse  increases,  more  radiation  is  scattered  into  the  higher  harmonics. 

Backscattered  Radiation.  In  the  backscattered  direction,  only  the  fundamental,  n  =  1, 
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is  nonzero.  In  the  limit  0—^0,  — >  1/2  and  Ji{ot)  — »  a/2.  Hence 


d?In  ^  e'^koNoM^o 
dudSl  0=0  4(1  +  a^/2)  ^ 

where  Gi{u)  is  given  by  Eq.  (41)  with  n  =  1. 


(50) 


IV.  RADIATION  PROPERTIES 


A.  Radiated  Power 

The  power  radiated  by  a  single  electron,  Pa,  undergoing  relativistic  quiver  motion  in 
an  intense  laser  field  can  be  calculated  from  the  relativistic  Larmor  formula  [24] 


*  3c  ^  \  dt  /  \dt  ) 


(51) 


Assuming  the  electron  orbit  is  a  function  of  only  the  variable  rj  =  z-\-  ct, 

2  /  j.  \  2l 


Ps  =  le^c('y  +  uaf 


Y 

)  \dv) 


(52) 


Using  the  orbits  described  in  Sec.  II,  the  power  radiated  by  an  electron  in  the  presence  of 
a  circularly  or  linearly  polarized  radiation  field  is  given  by 


Pa  ~  ^e^chlk^al  ■ 

(  sin^  koT],  linear, 


f  1/2. 


circular. 


(53) 


where  ho  is  given  by  Eq.  (12).  Averaging  the  above  expression  over  a  laser  period,  the 
ratio  of  the  radiated  power  to  the  incident  laser  power,  P,/Po>  cam  be  written  ais 


P,/P„  ~  16r|A5/3r5, 


(54) 


where  =  e^/rrieC^  is  the  claissical  electron  radius. 

The  total  power  radiated  by  a  laser  pulse  passing  through  a  uniform  distribution  of 
electrons  with  a  constant  density  no  is  given  by  Pt  =  NePa,  where  Ne  —  noLoai  is  the 
total  number  of  electrons  interacting  with  the  laiser  pulse  at  a  given  time,  Lq  =  cti  is 
the  laser  pulse  length  and  ol  is  the  effective  cross-section.  Assuming  a  Gaiissiam  laser 
pulse,  a  =  (aoro/ri,)exp(-r^/r|,),  where  rx,  is  the  laser  spot  size  and  tq  is  the  minimum 
spot  size,  the  effective  cross-section,  ai,,  can  be  foimd  by  letting  oq  — >  d  in  Eq.  (54)  and 
integrating  P*  over  r.  One  finds 


e-beam, 

plaisma. 


(55) 


where  /p  =  (1  +  ao/4)/(l  +  Oq/^)-  Eq.  (55),  the  top  expression  holds  in  the  short 
pulse  (electron  beam)  limit,  i.e.,  ho  =  7o(l  +  0o),  and  the  bottom  expression  holds  in  the 
long  pulse  (plasma)  limit,  i.e.,  ho  =  (1  +  ao/2)^/^.  Hence,  the  total  scattered  power  by  a 
uniform  electron  density  no  is  given  by 

Pt/Po  =  (87r/3)r2Lono/pho.  (56) 

As  example,  a  no  =  10^°  cm“^  plasma  interacting  with  a  1  ps  laser  pulse  with  ao  =  5  gives 
PtIPq  =  1-4  X  10“®.  The  ratio  of  the  total  scattered  energy  to  the  laser  pulse  energy  is 
approximately  PtL/PqLq,  where  L  is  the  total  length  over  which  the  laser  pulse  interacts 
with  the  electrons. 


B.  Resonance  Function 


Several  properties  of  the  radiation  spectra  can  be  etscertained  by  examining  the  res¬ 
onance  function,  R{k,nko),  given  by  Eq.  (34).  The  function  R{k,nko)  is  sheirply  peaked 
about  the  resonant  harmonic  frequencies,  ujn,  defined  by  ife  =  0,  which  can  be  written  as 


fiMoUfo 

[l-l-MoA(l-cos0)i’ 


(57) 


where  n  is  the  harmonic  niunber  and  Mq  is  the  relativistic  doppler  upshift  factor.  For  a 
plasma,  /3i  =  0  and  Mq  =  1,  which  gives  u»„  =  mvo,  independent  of  6.  For  a  relativistic 
electron  beam  with  Mq  ^  1,  the  radiation  is  primarily  backscattered  into  small  angles, 
9^  -C  1.  Hence  Un  —  nMou;o/(l  +  Mo^^/4),  which  indicates  a  maximmn  frequency  in  the 
backscattered  direction  along  the  axis,  9  =  0.  The  change  in  frequency  Aw  with  respect 
to  a  change  in  angle  A9  is  given  by 

lAu^l  ^  \Mo{9A9  +  My2)\ 

Wn  “  (2  +  Mo0V2)  ’  ^  ^ 


assuming  Mq  »  1.  Alternatively,  Eq.  (58)  can  be  solved  to  give  the  angular  spread  A9 
about  9  over  which  a  given  bandwidth  Aw  about  Wn  may  occupy.  For  a  relativistic  electron 
beam  with  Mq  >  1,  two  angles  are  of  particular  interest.  It  is  shown  below  that  for  a 
linearly  polarized  laser,  the  radiation  intensity  for  the  higher  harmonics,  n  »  1,  is  centered 
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•  1 
about  6  =  0,  whereas  for  circular  polarization,  the  intensity  is  centered  about  9q  =  2/Mo'*. 

For  these  two  angles,  Eq.  (58)  implies 


A0  ~  ( 


(  for  6  =  0, 


70 


[  (Aw/wn),  for  6  =  6o, 


(59) 


where  Mo  ^  47o  /7^  has  been  used. 

The  intrinsic  (i.e.,  associated  with  the  radiation  from  a  single  electron)  frequency  width 
Aa;„  of  the  radiation  about  a  resonant  frequency  a;„  can  be  found  by  letting 
and  integrating  the  function  R(k,  nko)  over  6a;,  which  gives 


,oo 

Acjn  =  /  d(6u})R{k,  ko)  =  Un/nNo. 

J  — OO 


(60) 


Hence,  Aa;„/a;„  =  I/tiNq,  where  Nq  =  Lq/Xq  is  the  number  of  wavelengths  in  the  laser 
pulse.  Furthermore,  R{k,  nko)  Aa;n6(a;  —  a;„)  as  iVo  — »  oo.  The  angular  width  A0„ 
within  which  can  be  foimd  radiation  with  frequencies  in  Aa;„  about  Un,  for  a  single  har¬ 
monic  n,  is  given  by  inserting  Eq.  (60)  in  Eq.  (59), 


A0n  ^  < 

70 


(  (1/nNo)^/^  for  6  =  0, 


(61) 


for  6  =  6o- 


I  (1/niVo), 

Alternatively,  similar  expressions  can  be  obtained  by  letting  6  =  6'  +  66  and  integrating 
R  ^o]  over  66.  It  should  be  pointed  out  that  Eqs.  (59)  and  (61)  apply  to  relativistic 

electron  beams  with  Mo  ^  1.  For  plasmas,  the  angrilar  width  occupied  by  a  given  Aa; 
about  Un  must  be  determined  by  considering  the  full  functioned  form  of  the  radiation 
spectrxim,  Eqs.  (36)  and  (46),  not  just  the  resonance  function  R{k,nko). 


C.  Ultra-Intense  Behavior 


For  values  of  oq  -C  1,  the  scattered  radiation  will  be  narrowly  peaked  about  the 
fundamental  resonant  frequency,  u)\  =  a;o/ [1  — /3i(l -I- cos0)].  As  oq  approaches  unity, 
scattered  radiation  will  appear  at  harmonics  of  the  resonant  frequency  as  well,  u;„  =  na;i. 
When  ao  »  1,  high  harmonic  (n  1)  radiation  is  generated  and  the  resulting  synchrotron 
radiation  spectrum  consists  of  many  closely  spaced  harmonics.  Finite  electron  energy 


22 


spread  effects  can  broaden  the  linewidth  causing  the  radiation  from  the  various  harmonics 
to  overlap.  For  example,  a  finite  thermal  axial  velocity  spread  will  lead  to  overlap  when 
{Au/u}n)th  ~  f/iT',  where  {Au}/u}n)th  is  given  below  by  Eq.  (77).  Hence,  in  the  ultra-intense 
limit,  i.e.,  cq  1,  the  gross  spectriim  appears  broadband,  and  a  continuiun  of  radiation 
is  generated  which  extends  out  to  a  criticad  frequency,  Wc,  beyond  which  the  radiation 
intensity  diminishes.  The  critical  frequency  can  be  written  as  tJc  =  ncU}R,  where  tie  is 
the  critical  harmonic  number.  It  is  possible  to  calculate  ric  by  examining  the  radiation 
spectrum,  Eqs.  (36)  and  (46),  in  the  ultra-intense  limit,  oq  ^  1- 

Asymptotic  properties  of  the  radiation  spectrum  for  large  harmonic  numbers,  n  »  1, 
can  be  analyzed  using  the  relationships  [25] 

fV2 


Jn{nz) 
J'„{nz)  ~ 


-(1-2^) 
TT 


(62) 


irz 


(1  -  z^)^/'^K2/3{nx), 


where  \z\  <  1  and  is  a  function  of  no  and  6, 


X  =  In 


1  4-  (1  -  -  Ini  -  (1  - 


(63) 


and  Ki/3,  K2/3  are  modified  Bessel  functions.  In  particular,  for  nx  »  1, 


Ki/3  ^  K2/3  ^  (7r/2nx)exp(-nx), 


(64) 


and,  hence,  only  harmonic  radiation  with  nx  <,  1  will  contribute  significantly  to  the  spec¬ 
trum.  The  critical  harmonic  number  is  defined  as  ricXmin  =  1>  i  e.,  Uc  =  l/imtu)  where 
Xmin  is  the  minimum  value  of  Eq.  (63).  Furthermore,  dx/dz  <  0  and  the  minimum  of  x 
occurs  at  imax-  Typically,  for  Cg  »  1,  1  —  z^ax  ^  ^  (®3)  can  be  expanded  to 

yield,  to  leading  order,  x^.n  -  (1/3)  (l  -  z^ax)  •  critical  harmonic  number  is  given 
by  the  inverse  of  this  expression. 

1.  Circular  Polarization 


For  a  circularly  polarized  incident  laser  field,  z  =  a/n,  where  a  is  given  by  Eq.  (47), 


i.e.. 


{ao/\^)  sing 
ho  [1  -  /?i(l  -I-  cos 0)1  ■ 


(65) 


23 


For  a  fixed  value  of  cq  ^  1,  the  maiximuin  value  of  i  is  given  by  Zmai  =  (ao/\/2)/(l  4- 
ao/2)^/^,  and  occurs  at  an  angle  Oq  given  by 

cos 00  =  (M)  -  1)/(M)  +  !)•  (66) 

Inserting  this  value  of  Zmax  into  Eq.  (63)  gives,  for  eg  »  1,  Xmin  —  2\/2/3a^  and,  hence, 

Tic  -  Zalf2y/2.  (67) 

Furthermore,  radiation  at  the  harmonic  ric  will  be  scattered  in  the  direction  6  =  Bq,  where 
00  is  given  by  Eq.  (66).  The  frequency  of  the  radiation  scattered  in  the  direction  0  =  0o  is 
given  by 

u,(e  =  0o)  =  nu;o(Mo  +  l)/2.  (68) 

For  a  plasma,  Mq  =  1  and  0o  =  ±7r/2,  i.e.,  the  high  harmonic  radiation  will  be  scattered 
perpendicular  to  the  incident  laser  field.  For  a  relativistic  electron  beam  with  Mo  ^  1, 
00  —  2/Mq'  and  the  high  harmonic  radiation  is  nearly  backscattered.  Physically,  0o  is 
related  to  the  pitch  angle  of  the  electron  orbit,  |ux|/|uj|  ~  2y/2/Mo^^  ^  Qo/70)  assuming 
Cq  »  1  and  Mq  »  1. 

The  asymptotic  properties  (n  »  1)  of  the  radiation  spectra  can  be  readily  obtained 
from  Eqs.  (46)  and  (62).  In  the  idtra-relativistic  limit,  Og  1,  the  radiation  is  confined  to 
small  angles  66  about  the  optimum  angle  0o.  i.e.,  0  =  0o  +  66,  where  66^  1.  Assuming 

n  :»  1,  Cq  1  and  66^  1,  Eqs.  (46)  and  (62)  give 

(kjdn  ~^°ir^c{l  +  'yW)  [(1  + 72002)^^/3(0 +  ^i^2/3(0j  .  (6  ) 

where 


UJc 

(70a) 

(Mq  +  1) 

fjjc-ric  2 

(706) 

ao(Mo  4- 1) 

^  “  2(2Mo)i/2  ■ 

(70c) 

Equation  (69)  holds  for  arbitrary  values  of  Mo,  i.e.,  electron  beams  of  arbitrary  energies  as 
well  as  stationary  plasmas.  In  Eq.  (70a),  ric  =  3ao/2\/2  and  the  factor  (Mg  4- 1)/2  is  the 
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relativistic  doppler  upshift  for  radiation  scattered  at  the  optimum  angle  6q,  as  indicated 
by  Eq.  (68).  The  expression  for  7  follows  from  Eq.  (8b)  assuming  Cq  >  1.  In  deriving 
Eq.  (69),  Eq.  (62)  was  used  and  the  summation  was  approximated  by  an  integral,  i.e., 
R{k,  nko)  ~  1/iVo  and,  hence,  nx 

Notice  in  the  limit  60  =  0,  tPl/dudO.  ~  where  ^  =  tj/wc.  A  plot  of  the 

function  y(^)  =  is  shown  in  Fig.  7.  The  function  Y{^)  is  maximum  at  ^  =  1/2 

and  decreases  rapidly  for  ^  >  1.  Half  the  total  power  is  radiated  at  frequencies  u  <  a;c/2 
and  half  at  a;  >  Wc/S.  This  can  be  shown  by  integrating  (Pl/dud^l  over  frequency  and 
angle  [10],  i.e.,  integrating  the  expression  given  below  by  Eq.  (71b)  over  frequency. 

Equation  (69)  is  Nq  times  the  standard  result  [24]  for  the  synchrotron  radiation  spec¬ 
trum  emitted  from  an  electron  moving  in  an  instantaneously  circular  orbit  in  the  ultra- 
relativistic  limit  with  a  radius  of  curvature  p  =  Sy^c/uJc-  Severed  well-known  properties 
[24]  follow  from  Eq.  (69),  for  example 

"  ~  [1  ^  ,  ^na) 

dSl  ~  48c  (1  +  7W)V2  [  7(1+  '■ 

=  2^/Z-Noy-  r  (TW) 

C  Uc  J2U./U.. 


The  peak  intensity  is  of  the  order  Noe^'y/c  and  the  toted  radiated  energy  is  of  the  order 
Noe^'ya>cfc.  The  peedc  intensity  occurs  at  the  optimum  emgle  Oq,  i.e.,  69  =  0,  at  approxi¬ 
mately  the  critical  frequency,  u;  ~  Uc,  i.e.,  n  Oi  ric  =  Za\l2y/2.  For  heu-monics  below  Uc, 
(a;  -C  Wc),  the  radiation  intensity  increases  as  (a>/u;o)^^^,  and  above  Uc  (u;  u>c),  the 

radiation  intensity  decreases  exponentially,  i.e.. 


d^I 

du6^l  6$=o 


JV„^  [r{2/3)]S“ 

7r*C 


d^I  ..  3e2  2/a;\  {  2a;\ 

1=;  j  [-zrj 


u  u;c, 

U)  »  LJc- 


(72a) 

(726) 


Furthermore,  for  u;  «:  Wc,  the  scattered  radiation  at  a  fixed  frequency  is  confined  to  an 
angular  spread  A60  =  about  9o,  whereas  for  u  >  Uc,  A60  =  (u;c/3u;)^/2/7- 

The  average  angular  spread  for  the  frequency  integrated  spectrum  is  {66^)^/^  ~  I/7. 

As  an  example,  the  peak  intensity  in  the  transverse  direction  (0  =  ?r/2)  of  each 


harmonic,  u  =  nuo,  is  shown  in  Fig.  8  for  the  case  of  a  high  intensity,  circularly  polarized 
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laser  pulse  encountering  an  electron  in  a  dense  plzisma.  Plots  for  two  different  intensities 
are  shown,  oq  =  4  and  cq  =  6.  The  arrows  indicate  the  approximate  critical  harmonic 
number,  ric  —  Oq,  for  each  case.  Asymptotically,  ao  >  1,  this  curve  approaches  the  form 
y(0  =  shown  in  Fig.  7. 

2.  Linear  Polarization 

For  a  linearly  polarized  incident  laser  field  in  the  limit  ao  -C  1,  upshifted  radiation 
at  the  fundamental  frequency  is  generated  in  a  narrow  cone  about  the  backscattered  di¬ 
rection,  n  ~  2ir6l,  where  6c  ~  1/ho-  However,  in  the  limit  ao  1,  a  near-continuum  of 
high  harmonic  radiation  is  generated  and  the  emission  cone  about  backscattered  direction 
widens  [10].  In  particular,  in  the  vertical  direction,  (j)  =  tt/ 2  (the  direction  normal  to 
the  x-z  plane  which  contains  the  electron  orbit),  emission  is  confined  to  the  vertical  angle 
~  1/ho-  In  the  horizontal  direction,  <^  =  0  (in  the  plane  of  the  electron  orbit),  the  emis¬ 
sion  angle  widens  and  is  confined  the  horizontal  angle  dh  ~  oo/ho,  which  is  determined 
by  the  deflection  angle  of  the  electron  in  the  x~z  plane  [10] .  The  asymptotic  properties 
of  the  radiation  spectrum  can  be  analyzed  using  Eqs.  (36)  and  (62).  Letting  $  represent 
the  observation  angle  in  the  vertical  direction,  i.e.,  =  ir/2,  then  in  the  limits  ao  1 

and  n  »  1,  <  1  and  the  coefficients  C*  and  C*  occurring  in  Eq.  (36)  are  given  by 

Cj  Ci  and  ~  (ao/ho)^J/^(£z),  where  additional  terms  of  order  1/ao  have  been 

neglected  and  n  =  2£  -I- 1  »  1.  Here,  for  linear  polarization. 


The  asymptotic  spectrum  near  the  axis  can  be  found  by  using  the  asymptotic  properties  of 
the  Bessel  functions,  Eq.  (62).  Notice  that  for  =  0,  Xjnax  —  8/3aQ.  Hence,  Ic  —  l/x„,ai, 
and  the  critical  harmonic  number,  Uc  2lc  is  given  by 

Uc  3a^/4.  (74) 


Using  Eqs.  (36)  and  (62),  the  asymptotic  spectrum  is  given  by 


<PI  12c2 
No 


dutdCl 


TT^C  (1  +  7^02) 


(75) 
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where 


(76a) 

(766) 

(76c) 


C=ii(l  +  ^79Y'^ 

LOc 

Uc  =  UcMoUq, 

7  =  ^o/2. 


In  deriving  Eq.  (75),  ^„i2(A:,  nfco)  l/iVo  and  lx  — >  C-  Several  subsequent  properties  of 
the  asymptotic  spectrum  follow  from  Eq.  (75).  As  was  the  case  for  circular  polarization, 
Eqs.  (71)- (72)  apply,  with  Nq  — >  4iVo,  66  6,  'y  j  and  where  Uc  is  given  by  Eq.  (76b). 

In  particular,  radiation  with  a;  ~  tJc  is  confined  to  a  vertical  angle  0^  —  1/7-  In  the 
horizontal  direction,  emission  is  confined  to  the  angle  ~  ao/j,  i.e.,  6h  ~  ao/7o  for  an 
electron  beam  and  ~  7r/2  for  a  plasma. 

As  an  example,  the  peak  intensity  on  axis  {6  =  0)  of  the  odd  harmonics,  u  =  tiMoUq,  is 
shown  in  Fig.  9  for  the  case  of  a  high  intensity,  linearly  polarized  laser  pulse  encountering 
a  counterstreaming  relativistic  electron  (70  =  5).  Plots  for  two  different  intensities  are 
shown,  ao  =  4  and  oq  =  6.  The  arrows  indicate  the  approximate  critical  harmonic  ntunber, 
Tic  —  3ao/4,  for  each  case.  Note  that  the  harmonic  intensity  is  plotted  versus  the  normalized 
frequency  uj/4'yQUo  ~  1.5ao.  Asymptotically,  oq  »  1,  this  curve  approaches  the  form 
r(0  =  f=jr|/3(e),  shown  in  Fig.  7. 
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V.  NON-IDEAL  EFFECTS 
A.  Electron  Energy  Spread 

The  above  analysis  has  assumed  ideal  electron  distributions,  i.e.,  thermal  and  en¬ 
ergy  spread  effects  have  been  neglected.  These  effects  are  important  in  determining  the 
frequency  line  width  of  the  scattered  radiation  [10].  For  example,  the  resonance  func¬ 
tion  R{k,nko)  indicates  that  if  a  thermal  axial  velocity  spread  Avth  is  introduced,  i.e., 
Pz  =  00  +  where  A0th  =  ^vth/ct  then  the  scattered  radiation  along  the  axis  will  be 
shifted  in  frequency  away  from  by  Autk,  where 

(AuMth  =  27^A0th-  (77) 

For  a  plasma,  A0th  is  related  to  the  initial  plasma  thermal  energy,  Eth,  by  A0th  = 
{2Eth/TneC^y^^.  For  an  electron  beam,  A0th  is  related  to  the  initial  normalized  energy 
spread,  A7/70,  by  A0th  =  A^IjqPq.  As  an  example,  a  plasma  with  a  temperature  of  100 
eV  would  produced  a  thermal  bandwidth  of  {Au/un)ih  —  4%. 

In  actual  electron  beams,  the  electrons  may  have  an  average  angular  spread  as  well  as 
an  average  energy  spread,  represented  by  emittance  and  intrinsic  energy  spread,  respec¬ 
tively.  The  normaUzed  beam  emittance  is  given  by  en  =  7o»’6^6>  where  is  the  average 
electron  beam  radius  and  6b  is  the  average  electron  angtilzir  spread.  The  fractional  lon¬ 
gitudinal  beam  energy  spread  due  to  emittance  is  {AE/Eb)t  =  e^/2r^,  where  Eb  is  the 
initial  beam  energy.  Electron  beams  may  also  have  an  intrinsic  energy  spread,  {AE/Eb)i, 
due  to  various  reaisons,  such  as,  voltage  variation,  finite  pulse  length  effects,  etc.  The  total 
spectral  width  of  the  radiation  about  the  harmonic  Un  is 

{Auj/ujnh  ^  [(Au;/u;n)^  -I-  (Atj/wn)*  +  (Au;/u;„)?]  ,  (78) 

where  (Au;/u;„)o  =  I/tiNq  is  the  finite  interaction  length  spectral  width  contribution, 
{Au}/ii}n)t  =  Cn/^b  the  emittance  broadened  spectral  width  and  (Aa;/u;„)i  =  2{AE/Eb)i 
is  the  intrinsic  energy  spread  broadening  contribution.  The  radiation  with  total  spectral 
width  {Au}lu}n)T  is  confined  to  the  angle  Bt  —  (Au;/u;„)^^^/7o.  This  consequently  reduces 
the  spectral  intensity,  (Pi /dudU,  of  the  scattered  radiation  from  an  electron  beam  for  a 
particular  harmonic  by  approximately  B\IB\. 
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If  a  particular  application  requires  a  bandwidth  {Au)/iJn)s  <  1,  this  radiation  can  be 
foTind  within  the  angle  0^,  where 

=  [(^‘*^/‘*^n)s  +  (^^/^n)T)]  /lo- 

If  a  bandwidth  {Au/un)s  ^  (^tj/wn)r  is  required,  all  the  radiation  within  a  cone  of 
half-angle  6^  c:^  6s  =  {Au/unis^ /'^o  used.  To  obtain  a  bandwidth  (Au/uJn)s  ^ 

{A(jj / Un)T  1  the  radiation  within  the  cone  ^  dx  =  (Aa;/a;„)y^^/7o  must  be  filtered  using 
a  monochromator.  As  an  illustration,  for  an  rf  linac  electron  beam  with  e„  ~  5  mm-mrad, 
rj,  =  50  ^m  and  70  =  100,  {AE/Eb)t  ~  0.5%  and  {Au/u}n)c  —  1%-  Since  the  intrinsic 
energy  spread  is  typically  ~  1%  and  Nq  >  300,  the  total  spectrail  width  of  the  unfiltered 
LSS  radiation  is  typically  {Au/uJn)T  —  1%  and  is  confined  to  the  angle  6t  —  1  mrad. 

B.  Electron  Beam  Energy  Loss 

As  the  electron  beam  radiates  via  nonlinear  Thomson  scattering,  the  electron  beam 
will  lose  energy.  The  rate  of  loss  of  electron  beam  energy  is  equad  to  the  scattered  power, 
meC^d'y/dt  =  —Pa,  where  P*  is  given  by  Eq.  (54).  Assuming  ho  ~  47^,  the  electron  beam 
energy  will  evolve  [5]  according  to  7  =  7o/(l  +  tfTR),  where  t  is  the  electron  beam-laser 
interaction  time  and  tr  =  3/(4crefcoao7o)»  where  a  linearly  polarized  laser  field  has  been 
assumed.  In  practical  units,  this  can  be  written  as 

tr[ps]  ~  1.6  X  10^^ [MeV]  ^W/cm^j  .  (80) 

One  consequence  of  the  loss  of  electron  beam  energy  is  the  introduction  of  an  additional 
source  of  enhanced  bandwidth,  (Au;/a;„)R  =  2(70  —  7)/70)  where  70  —  7  =  jt/rR.  For 
typicad  values  of  laaer  p\ilse  lengths  and  intensities  of  interest,  I/tr  «:  1,  and  this  effect 
is  smadl.  As  an  example,  a  2  ps  (t  =  1  ps)  laser  pulse  with  intensity  Iq  =  2.6  x  10^^ 
W/cm^  (ao  =  0.43)  interacting  with  a  Ej,  =  40  MeV  (70  =  79)  electron  beam  gives 
{Au)/u}n)R  —  0.13%. 

C.  Ponderomotive  Density  Depletion 

In  a  high  density  plasma,  the  tramsverse  ponderomotive  force  from  the  radial  gradients 
in  the  laser  pulse  profile  cam  displace  the  plasma  electrons  leaiding  to  a  density  depression 
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on  axis.  In  the  long  pulse  limit,  the  density  depression  can  be  calculated  by  equating  the 
electrostatic  force  with  the  ponderomotive  force,  =  Vj.7x,  which  is  the  adiabatic 
response  of  the  plaisma  electrons  to  the  transverse  ponderomotive  force  [23,26].  This  gives 
an  equilibrium  density  profile  of 

Tie/no  =  1  +  fcp  ^V^(l  -t-  a^/2)^/^,  (81) 


where  Ue/no  >  0  has  been  assumed.  Assuming  a  Gaussian  transverse  profile  of  the  form 
|a|  ~  exp(— r^/r^),  Eq.  (81)  indicates  that  the  density  along  the  axis  is  given  by 


ne(r  =  0)  ^  aQ>i 
riQ  2iT^rl 


(82) 


where  Ap  =  27r/fcp.  As  an  example,  a  high-density  plasma  with  ro  =  15  ^m,  Ap  =  5  /xm 
and  oo  =  7  gives  a  density  depression  along  the  axis  of  And  no  =  5%.  This  density 
depression  reduces  the  total  number  of  electrons  scattering  radiation,  hence,  the  total 
scattered  power  ~  will  be  reduced.  Furthermore,  in  a  high  density  plasma,  the  ef¬ 
fects  of  relativistic  self-focusing,  which  occurs  for  pump  laser  powers  above  a  critical  power, 
Pc[GW]  ~  17(Ap/Ao)^,  along  with  the  effects  of  a  density  depletion  on  axis,  can  provide  op¬ 
tical  guiding  and  significantly  extend  the  laser-plasma  interaction  distance  [5,22,23,26,27]. 
For  a  relativistic  electron  beam  in  the  short  pulse  limit,  tl  -C  7o/t^p>  the  magnitude 
of  the  electron  density  perturbation,  Aue,  due  to  the  ponderomotive  force,  is  given  by 
[Ane/noj  {Loao/'yoro)^  -C  1,  consistent  with  the  discussion  at  the  end  of  Section  II. 


D.  Plasma  Dispersion 


The  frequency  of  the  scattered  radiation  cam  be  affected  by  the  dispersion  properties 
of  electromagnetic  radiation  in  a  plasma.  In  the  long  pulse  Umit,  the  nonhnear  dispersion 
relation  for  radiation  of  frequency  u  and  wavenumber  k  is  given  [5]  by  ~  -t-u;p/7j.. 

Notice  that  the  dispersion  relation  is  different  for  radiation  within  the  region  of  the  pump 
laser  pulse,  7x  =  (1  +  ao/2)^/^,  and  for  radiation  propagating  in  the  plasma  outside  of  the 
pump  laser  pulse,  7j.  =  1.  In  particular,  for  backscattered  radiation,  the  radiation  will 
transit  a  counterstreaming  boundary  region  at  the  trailing  edge  of  the  pump  laser  pulse. 
As  the  backscattered  radiation  transits  this  boundary  region,  coimterstreaming  at  the 
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group  velocity  of  the  pump  laser  pulse,  Vg,  the  frequency  and  wavenumber  of  the  scattered 
radiation  will  be  shifted  [28].  Hence,  the  detected  frequency,  Ud,  of  the  backscattered 
radiation  will  be  shifted  from  the  frequency  at  which  it  is  scattered,  u;,  within  the  laser 
pulse.  The  detected  frequency  uid  is  related  to  the  scattered  frequency  u  by  requiring  the 
phase  of  the  scattered  radiation  to  be  continuous  across  the  boundzu'y  at  the  trailing  edge 
of  the  laser  pulse  [28],  u  +  Vgk  =  (jjd  +  Vgkd,  where  Vg  =  c(l  —  and  a  square 

laser  pulse  profile  has  been  assumed  for  simplicity.  Using  the  dispersion  relation  to  solve 
for  k  and  kd  in  terms  of  u  and  oJd,  respectively,  and  assuming  -C  1,  implies 


U) 


(83) 


Hence,  for  backscattered  radiation,  the  detected  frequency  will  be  upshifted  from  the  scat¬ 
tered  frequency.  Furthermore,  depletion  of  the  electron  plasma  density  within  the  region 
of  the  laser  pulse  by  the  transverse  ponderomotive  force  will  produce  a  additional  upshift 
for  similar  reasons.  This  effect  can  be  approximated  by  replacing  l/‘y±  with  ne/7x”o 
in  Eq.  (83),  where  Ue/no  is  given  by  Eq.  (82).  The  maximum  frequency  upshift  for  the 
backscattered  radiation  can  be  estimated  by  Au>a/w  Ci;  which  is  typically  small. 

Radiation  scattered  in  the  transverse  or  forward  directions  will  not  experience  a  frequency 
shift  by  these  mechanisms. 
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VI.  LASER  SYNCHROTRON  SOURCES 
Nonlinear  Thomson  scattering  can  be  used  as  a  mechanism  for  generating  x-ray  radi¬ 
ation  [1-9].  In  such  a  laser  synchrotron  source  (LSS),  intense  laser  pulses  are  backscattered 
from  a  counterstreaming  relativistic  electron  beam  or  from  a  dense  plasma  [5-7].  The  LSS 
has  the  potential  for  providing  a  compact  source  of  tunable,  short  pulse  radiation,  in  the 
soft  to  hard  x-ray  regime.  Two  examples  of  LSS  configurations  will  be  discussed,  one 
using  a  relativistic  electron  beam  to  generate  hard  x-rays  (30  keV,  0.4  A),  and  the  other 
using  a  dense  plasma  to  generate  soft  x-rays  (300  eV,  40  A).  In  the  electron  beam  LSS, 
short  wavelengths  are  generated  by  exploiting  relativistic  doppler  factor,  i.e.,  A  =  Ao/47o, 
assuming  7o  2>  1  and  ag  <C  1.  In  the  plasma  Ll  ^  short  wavelengths  are  generated  by 
exploiting  the  nonlinear  harmonic  factor,  i.t  a  =  Xq/uc,  where  Uc  ~  Oq  2>  1  is  assumed. 
Both  configurations  will  utilize  the  recently  developed  sohd-state  laser  technology  based 
on  chirped-pulse  amplification  (CPA)  [18-20].  Lasers  based  on  CPA  are  relatively  com¬ 
pact  systems  capable  of  dehvering  ultrahigh  powers  (  ^  10  TW)  and  intensities  (  >  10^® 
W/cm^)  in  ultrashort  pulses  ( ;$  1  ps).  Currently,  the  repetition  rates  of  TW  CPA  systems 
are  limited  to  <  10  Hz  [19,20].  A  summary  of  the  current  state-of-the-art  in  CPA  laser 
technology  can  be  found  in  Ref.  [20]. 

A.  Electron— Beam  LSS 


An  electron-beam  LSS  configuration  consists  of  backscattering  a  linearly  polarized 
laser  pulse  from  a  counterstreaming  relativistic  electron  beam.  Two  important  quantities 
cheu’acterizing  the  resulting  synchrotron  radiation  are  the  photon  flux,  F,  defined  as  the 
number  of  photons  per  second  within  a  specified  bandwidth,  and  the  photon  brightness, 
B,  defined  as  the  phase  space  density  of  the  photon  flux.  The  intensity  distribution  for 
backscattered,  0  =  0,  radiation  at  the  fundamented,  n—  1,  in  the  limit  Cq  C  1  and  7o  >  1 
(i.e.,  u;  ~  (D  =  47oa;o),  is  given  by 


d^7(0) 


Sttc^  ~ 


sin(7r(u;  —  u})No/ij) 
ir{uj  —  {I))Nq/G} 


(84) 


ais  indicated  by  Eq.  (39).  The  angular  density  of  the  flux,  dF/dil,  i.e.,  the  peak  number  of 
photons  in  a  specified  frequency  range  ui  <  u  <  uJ2  emitted  per  second  per  unit  sohd  angle 
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by  the  micropulse  in  the  forward  direction,  can  be  determined  from  Eq.  (84)  by  integrating 
over  the  frequency  range  Aus  =  UJ1—U2,  multiplying  by  the  electron  flux  interacting  with 
the  laser,  and  by  dividing  by  the  energy  per  photon,  hu>.  The  electron  flux  interacting 
with  the  laser  field  is  given  by  Nb  =  fh/e,  where  Ib  is  the  peak  micropulse  current  and  / 
is  the  filling  factor,  i.e.,  /  =  ffo/o’6  for  cq  <  Cb  and  /  =  1  for  (Tq  >  (Tb,  where  oq,  ab  are 
the  cross-section  areas  of  the  laser  and  electron  beam,  respectively.  The  angular  density 
of  the  flux  is  given  by 


^  =  aN^N^alil  ■  \ 


{  Nq{AuIq)s,  for  (Alj/u>)s  «:  I/Nq, 


(85) 


1,  for  (Aa;/a))s  »  I/Nq, 

where  a  =  1/137  and  Fq  denotes  the  spectral  flux  for  an  ideal  electron  beam,  i.e.,  zero 
emittance  and  energy  spread.  For  an  ideal  electron  beam,  the  spectral  flux  with  spectral 
width  {Au/u})s  is  given  by  Fq  ~  27rfl^(dFo/dn),  where  =  6l  +  i.e.. 


Fo  cii  2TroiNoNb0.l{Au/Q)s, 


(86) 


which  is  valid  for  all  values  of  (Aa)/u;)s  <  1.  For  a  realistic  electron  beam  with  finite 
emittance  and  energy  spread,  the  photon  flux,  F,  is  identical  to  the  the  ideal  case,  i.e., 
F  =  Fo.  The  angular  density  of  the  flux,  dF/dCl,  however,  is  reduced,  since  the  photons 
are  now  spread  out  over  a  larger  radiation  angle  0s >  where  0s  is  given  by  Eq.  (79),  i.e., 
dF/dn  ~  Fo/27r0|:. 


The  spectral  brightness  is  the  phase  space  density  of  F.  Hence,  B  =  F/(27r)^(/l0s)^, 
where  (i20s)^  is  the  phase  space  area  of  the  photon  beam.  The  quantity  R  is  the  total 
effective  size  of  the  radiation  source  and  is  given  by  R?  ^  +  {0EiLf47r)^,  where  9^,  = 

6\  +  0?,  6i  =  {Au)/u))y^ /'fo,  and  r,  is  the  smaller  of  r*,  and  ro/2.  Here  L  is  the  laser- 
electron  interaction  distance.  The  spectral  flux  and  brightness  for  a  non-ideed  electron 
beam,  in  terms  of  practical  units,  are  given  by 
photons^ 


sec 


=  8.4  X  10^^f{L/ZR)Ib[A]Po[GW]{Au}/u>)s, 


(87a) 


photons 


sec  •  mm^  •  mrad"^ 


=  8.1  X  10®/(L/ZR)(/6[A]/r^[mm]) 

(Au>/u))5/(1  +  6) 


F6^[MeV]Fo{GW] 


{Au}/u})s  +  {Au/iI))T 


(876) 
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where  6  =  (0E<L/47rrs)^  is  typically  1.  The  interaction  length  is  the  smaller  of  twice  the 
Rayleigh  length  {Zr  =  ■ktI/Xq)  or  one-half  the  laser  pulse  length,  i.e.,  L  =  min  [2Zr,  Lo/2], 
unless  it  is  further  hmited  by  the  specific  geometry  of  the  experiment. 

As  an  example,  consider  an  electron  beam  LSS  which  generates  0.4  A  (30  keV)  x- 
rays.  For  a  Aq  =  1  /xm  incident  laser,  A  =  Ao/47o  =  0.4  A  imphes  70  =  79  {Eb  —  40 
MeV),  assuming  Ug  C  1.  A  CPA  laser  will  be  assumed  with  tq  =  2  ps,  Pq  =  10  TW  and 
To  =  50  fj,m,  which  imphes  Iq  =  2.6  x  10^^  W/cm^,  cq  =  0.43  and  Zr  =  7.9  mm.  An 
electron  beam  from  an  rf  linac  wiU  be  assumed  with  peak  current  Ib  =  200  A,  micropulse 
duration  Lb/c  =  1  ps,  beam  radius  rb  =  50  /im,  energy  spread  (AEfEb)  =  0.5%  and 
normaUzed  emittance  €„  =  5  mm-mrad.  The  interaction  length  is  one  half  the  laser  pulse 
length,  L  =  300  /2m,  and  the  x-ray  pulse  duration  is  the  micropulse  duration,  =  1  ps. 
The  effective  bandwidth  is  {Au}/u})t  —  1.4%  and  this  radiation  is  confined  to  a  cone  angle 
of  6t  —  1-5  mrad.  The  total  flux  with  (Au/cj)s  ~  1  within  the  cone  6c  ~  I/70  ~  12  mrad  is 
F  =  6.4  X  10^^  photons/sec.  The  peak  brightness  with  {Au}/Q)s  =  0.1%  is  5  =  2.9  x  10^® 
photons/s-mm^-mrad^.  The  parameters  for  this  electron  beam  LSS  are  summarized  in 
Table  I. 

For  simpUcity,  a  counterstreaming  laser-electron  beam  geometry  has  been  assumed  in 
which  the  x-ray  pulse  length  is  approximately  the  electron  micropulse  length.  Shorter  x- 
ray  pulse  lengths  can  be  obtained  by  either  reducing  the  laser  Rayleigh  length  or  changing 
the  laser-electron  beam  intersection  angle  [6,8].  (Kim  et  al.  [8]  have  suggested  scattering 
at  90°  to  obtain  ultrcishort  x-ray  pulses.)  In  principle,  both  these  methods  may  lead  to 
the  production  of  ultrashort  x-ray  pulses,  with  pulse  durations  on  the  order  of  the  laser 
pulse  duration. 


B.  Plasma  LSS 

To  produce  x-rays  with  a  Aq  ~  1  uva.  laser  and  a  stationary  plasma,  it  is  necessary 
to  use  xiltrahigh  intensities,  Cg  1.  Nonlinear  Thomson  scattering  will  then  occur  in  the 
asymptotic  limit,  in  which  a  neau:  continuum  is  produced  with  harmonics  extending  out 
to  the  critical  harmonic  number,  Uc  ~  Og,  as  discussed  in  Section  2.  Consider  a  Unear 
polarized  laser  with  Og  >  1  interacting  with  a  dense  plasma.  In  the  near  backscattered 
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direction  the  radiation  spectrum  scattered  by  a  single  electron  is  given  by 

iM  -  N 

duda  ^°2ir2c 


as  indicated  by  Eq.  (75),  where  Y  =  i  —  =  UcIjJq  and  Uc  —  3ao/4. 

For  a  collection  of  electrons  in  a  plasma,  the  totEil  energy  radiated  is  given  by  Et  = 
iVeJ(O),  where  =  rieffoLp  is  the  total  number  of  electrons  with  which  the  laser  interacts, 
rig  is  the  plasma  electron  density,  (Tq  =  7rr^/2  is  the  laser  cross-section  and  Lp  is  the 
laser-plasma  interaction  distance.  Typically,  Lp  ~  2Zr  =  27rro/Ao,  assuming  vacuum 
di&action.  The  effects  of  relativistic  optical  guiding,  however,  could  substantially  increase 
the  interaction  distance  [5,22,23,26,27].  Geometric  arguments  indicated  that  the  x-ray 
pulse  length  in  the  backscattered  direction  is  given  by  Lx  ~  2Lp(l  +  Ll/ALp)^^^  ~  2Lp, 
where  Lq  is  the  laser  pulse  length  and  Lq/4Lp  •C  1  has  been  assumed.  The  total  power  in 
the  backscattered  direction  is  Pt  =  cEtILx  and  the  photon  flux  is  P  =  PtIJuj}.  Hence, 
the  flux  intensity,  defined  to  be  dF/dCl,  for  photons  in  the  frequency  range  Aus  about  u 
in  the  near  backscattered  direction,  is  given  by 


dFJdO.  ^  (3ac/87r)JVongroOo(Au;/u;)5y(u>/u;c)-  (89) 


Recall  that  the  solid  angle  over  which  the  photons  with  frequencies  near  Uc  are  scattered  is 
relatively  large,  i.e.,  ~  2\/2/oo  in  the  vertical  direction  and  9h  ~  7r/2  in  the  horizontal 

direction.  The  total  photon  flux,  F,  can  be  estimated  by  multiplying  Eq.  (89)  by  the 
appropriate  solid  angle  over  which  the  photons  eu’e  to  be  collected.  The  brightness,  B,  of 
the  backscattered  photons  can  be  estimated  by  B  ~  {dFfdd)/nrQ.  In  practical  units,  the 
photon  flux  intensity  and  brightness  are  given  by 


dF  photons 
dn  s  •  mrad^ 


~  3.65  X  10  ^To  [ps]  Ao  [^imj  Ug  [cm  Po  [TW] 

■  {Au}/u))sY{ijjl(jJc), 


photons 
mm^  •  mrad^ 


~  1.80  X  10  [ps]  Ao[^m]ng  [cm  ^]  /q  [w/cm^j 

•  (Au;/u;)sy(u;/u>c), 


(90a) 

(905) 


As  an  example,  consider  a  plasma  LSS  which  generates  40  A  x-rays.  For  a  Ao  =  1 
To  =  1  ps  incident  laser  pulse,  A  =  Ao/ug  =  40  A  implies  rig  =  250  and  oq  =  6.9,  which 
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corresponds  to  a  laser  intensity  of  Jo  =  6.6  x  10^®  W/cm^.  Assuming  a  laser  spot  size 
of  To  =  15  /im  gives  a  laser  power  of  Pq  =  230  TW  and  a  laser-plasma  interaction  of 
length  of  Lp  ~  2Zr  =  1.4  mm.  The  x-ray  pulse  duration  is  Tx  ~  2Lp/c  =  9.4  ps.  A 
plasma  density  of  rie  —  10^°  cm~^  implies  a  flux  intensity  of  dF/dfl  ~  2.1  x  10^®(Aa;/a;)5 
photons/s- mrad^  and  a  brightness  of  S  =  2.9  x  10^^(Au;/u;)s  photons/s-mm^-mrad^.  The 
parameters  for  this  plasma  LSS  are  summarized  in  Table  II. 

For  simplicity,  the  generation  of  backscattered  {6  =  0)  x-rays  from  the  interaction  of 
a  linearly  polarized  laser  and  a  plasma  has  been  considered.  For  this  case,  the  x-ray  pulse 
length  is  of  the  order  of  a  few  Rayleigh  lengths.  However,  Eqs.  (36)  and  (46)  indicate  that 
somewhat  leirger  fluxes  of  x-rays  are  emitted  in  the  transverse  direction  (0  =  7r/2)  for  both 
circularly  and  linearly  polarized  lasers  incident  on  a  plasma.  Hence,  by  collimating  the 
transverse  emission  from  a  plasma,  ultrashort  x-ray  pulses  can  be  obtained  with  durations, 
in  principle,  on  the  order  of  the  laser  pulse  duration. 
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VII.  CONCLUSION 


A  comprehensive  theory  describing  the  nonlinear  Thomson  scattering  of  intense  laser 
fields  from  beams  and  plasmas  has  been  presented.  This  theory  is  vahd  for  hnearly  or  cir¬ 
cularly  polarized  incident  laser  fields  of  arbitrary  intensities  and  for  electrons  of  arbitrary 
energies.  Exphcit  expressions  for  the  intensity  distributions  of  the  scattered  radiation 
were  calculated  anadytically  and  evaluated  numerically.  The  space-charge  electrostatic 
potential,  which  is  important  in  high  density  plaismas  and  prevents  the  aodal  drift  of  elec¬ 
trons,  wais  included  self-consistently.  Various  properties  of  the  scattered  radiation  were 
examined,  including  the  linewidth,  angular  distribution,  and  the  behavior  of  the  radia¬ 
tion  spectra  at  ultrahigh  intensities  (cq  1).  Non-idecil  effects,  such  as  electron  energy 
spread  and  beam  emittance,  which  can  broaden  the  linewidth  and  angular  distribution  of 
the  scattered  radiation,  were  disciissed.  These  results  were  then  applied  to  possible  LSS 
configurations. 

The  general  formula  for  the  frequency  of  the  Thomson  backscattered  {6  =  0)  radiation 
is  given  by  Un  =  tiMoUq,  where  n  is  the  harmonic  number  and  Mq  is  the  doppler  multiph- 
cation  factor,  given  by  Eq.  (17).  For  a  linearly  polarized  laser,  only  odd  harmonics  exist 
in  the  backscattered  direction,  whereas  for  circular  polarization,  only  the  fundamental  is 
nonzero  in  the  backscattered  direction.  Both  odd  and  even  harmonics  can  exist  at  off-axis 
angles.  General  expressions  for  the  scattered  intensity  distributions  are  given  by  Eqs.  (36) 
and  (46).  Generation  of  x-rays  at  short  wavelengths  require  Mo  ^  1  and/or  n  1.  The 
intrinsic  linewidth  (i.e.,  for  a  cold  electron  distribution)  of  a  particular  harmonic  is  given 
by  Au/uji  =  l/nNo,  where  Nq  is  the  number  of  laser  periods  with  which  the  electrons 
interact.  Since  Nq  >  300,  smzdl  linewidths  can  be  achieved.  Non-ideal  effects,  such  as 
energy  spread  and  beam  emittance,  can  broaden  the  linewidth,  as  indicated  by  Eq.  (78). 
When  al  <SC  1,  radiation  is  scattered  only  at  the  fundamental.  When  Cq  2>  1,  a  multi¬ 
tude  of  hairmonics  are  produced,  which  results  in  a  near-continuum  of  scattered  radiation 
extending  out  to  a  critical  harmonic  number,  ric  ~  Cp,  beyond  which  the  intensity  of  the 
radiation  rapidly  diminishes.  Expressions  for  the  scattered  intensity  distributions  in  the 
ultra-intense  limit  are  given  by  Eqs.  (69)  and  (75).  The  polarization  of  the  scattered  radi¬ 
ation  can  be  adjusted  by  changing  the  polarization  of  the  incident  laser.  Scattering  from 
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an  electron  beam  has  the  additional  advantage  of  well-collimated  radiation.  For  70  ^  1 
and  Oq  -C  1,  the  upshifted  radiation  is  confined  to  a  cone  about  the  backscattered  direc¬ 
tion  of  half-angle  6  ~  (Aa;/d))^/^/7o.  Scattering  from  a  plasma  has  the  advantage  in  the 
attainability  of  high  electron  densities,  the  photon  flux  and  brightness  scaling  linearly  with 
density. 

A  LSS,  based  on  the  nonlinear  Thomson  scattering  of  intense  lasers  from  electron 
beams  or  plasmas,  may  provide  a  practical  method  for  producing  x-ray  radiation.  The 
LSS  has  a  number  of  potentially  unique  amd  attractive  features  which  may  serve  a  variety 
of  x-ray  spectroscopic  and  imagining  applications.  These  features  include  compactness, 
relatively  low  cost,  tunability,  narrow  bandwidth,  short  pulse  structure,  high  photon  energy 
operation,  well-collimated  photon  beams,  polarization  control,  and  high  levels  of  photon 
flux  and  brightness.  Specific  examples  of  an  electron-beam  LSS  and  a  plasma  LSS  were 
given,  as  svunmarized  in  Tables  I  and  II.  An  electron-beam  LSS,  designed  to  generate  30 
keV  (0.4  A)  photons  with  a  Aq  =  1  fim  laser  with  cq  <  1,  requires  a  40  MeV  electron 
beam  (approximately  300  times  lower  energy  electrons  than  required  by  a  conventioneil, 
storage-ring  synchrotron).  This  electron  beam  LSS  generates  1  ps  x-ray  pulses  with  a  high 
peak  flux  (  >  10^^  photons/s)  and  brightness  (  >  10^®  photons/s-mm^-mrad®,  0.1%BW).  A 
plasma  LSS,  designed  to  generated  40  A  (0.3  keV)  photons  with  a  Aq  =  1  fim  laser,  requires 
ao  =  6.9  {Jo  =  6.6  x  10^®  W/cm^).  This  plasma  LSS  generates  <  10  ps  x-ray  pulses  with 
a  high  peak  flux  (  >  10^^  photons/s,  10^  mrad^)  and  brightness  (  >  10^®  photons/s-mm^- 
mrad^,  0.1%BW).  These  peak  values  of  flux  and  brightness  compare  favorably  to  those 
obtained  in  conventional  synchrotrons.  High  levels  of  average  flux  emd  brightness  are 
presently  limited  by  laser  technology.  The  recent  advances  in  compact,  solid-state  lasers, 
based  on  chirped-pulse  amplification,  are  capable  of  generating  the  ultrahigh  intensities 
(oo  ^1)  needed  to  experimentally  explore  Thomson  scattering  and  LSS  x-ray  generation 
in  the  nonlinear  regime. 

This  paper  has  been  restricted  to  the  discussion  and  analysis  of  x-ray  generation  by 
the  Thomson  (incoherent)  scattering  of  intense  lasers  from  beams  and  plasmais.  However, 
for  sufl^ciently  cold  electron  distributions,  it  is  also  possible  to  generate  short-wavelength 
radiation  by  the  stimulated  (coherent)  backscattering  of  intense  lasers  from  beams  and  plas- 
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mas  [5,21,29,30].  Stimtilated  backscattered  harmonic  generation  may  provided  a  method 
for  producing  coherent  x-rays  via  a  laser-pumped  free  electron  laser  (LPFEL).  Advances 
in  CPA  lasers  and  in  high-brightness  electron  beams  may  soon  provide  the  necessary  tech¬ 
nology  to  realize  compact  sources  of  both  incoherent  (LSS)  and  coherent  (LPFEL)  x-rays. 
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Table  I 


Parameters  for  an  Electron  Beam  LSS 


Incident  Laser  Parameters 

Wavelength,  Aq 
Pulse  Length,  Lq/c 
Pe£ik  Power,  Pq 
Intensity,  Iq 
Strength  Parameter,  ao 
Spot  size,  ro 
Rayleigh  Length,  Zr 

Electron  Pulse  Parameters 

Beam  Energy,  Et 
Beam  Current,  /{, 

Beam  Pulse  Length,  Li/c 
Beam  Radius,  Tb 
Beam  Energy  Spread,  {AE/Eb)i 
Beam  Emittance,  Cn 

X-Rav  Pulse  Parameters 

Photon  Energy,  Ep 
Photon  Piilse  Length,  Lb/c 
Peak  Photon  Flux,®  F 
Photons/Pulse,®  FLbjc 
Peak  Brightness  (0.1%  BW),  B 
Angular  Spread,  0c  ~  1/7 


1  fim 

2  ps 
10  TW 

2.6  X  lO^"^  W/cm2 
0.43 
50  fim 
7.9  mm 


41  MeV 
200  A 
1  ps 
50  fxm 
0.5% 

5  mm-mrad 


30  keV 
1  ps 

6.4  X  10^^  photons/s 

6.4  X  10®  photons/pulse 

2.9  X  10^®  photons/s- mm^-mrad^ 

12  mrad 


®Includes  all  photons  within  the  ~  I/7  angle,  implying  ~  100%  BW. 
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Table  II 


Paramftters  for  a  Plasma  LSS 


Incident  Laser  Parameters 

Wavelength,  Aq 
Pulse  Duration,  Tq 
Peak  Power,  Pq 
Peak  Intensity,  Iq 
Strength  Parameter,  uq 
Spot  size,  ro 
Rayleigh  Length,  Zr 

Plasma  Parameters 

Electron  Density, 

Interaction  Length,  2Zr 

X-Rav  Piilsp  Parameters 

Wavelength,  A* 

Photon  Energy,  Ep 
Photon  Pulse  Length,  r* 

Flux  Intensity  (0.1%  BW),  dF/dQ 
Brightness  (0,1%  BW),  B 
Photon  Flux“  (100%  BW),  F 

“Includes  photons  with  (Aaj/a;)s  ~  1  withi 


1  fim 
1  ps 

230  TW 

6.6  X  10^®  W/cm2 

6.9 

15  urn. 

710  fxm 


102o  cm-3 

1.4  mm 


40  A 
310  eV 

9.4  ps 

2.1  X  10^®  photons/s-mrad^ 

2.9  X  10^®  photons/s-mm^-mrad^ 

6.5  X  10^^  photons/s 

a  solid  angle  dfl  ~  with  0  =  10  mrad. 
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Fig.  1  Schematic  diagram  showing  the  Thomson  scattering  of  an  intense  laser  field  from 
a  free  electron. 
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Fig.  3  The  normalized  intensity  at  the  harmonic  resonances,  u/wo  =  n,  as  &  function  of 
angle,  6,  in  the  <p  =  0  plane,  of  the  radiation  scattered  by  a  dense  plasma  electron 
from  a  linearly  polarized  laser  pulse  {No  =  7).  (a)  shows  the  first  three  harmonics 
for  ao  =  0.5,  (b)  shows  the  first  six  harmonics  for  oo  =  1.0,  and  (c)  shows  the 
first  twelve  harmonics  for  an  =  2.0. 
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Fig.  4  The  hanuonic  amplitude  function,  F„(oo),  as  a  function  of  (ao/4)/(l  +  Oo/2),  for 
the  first  ten  odd  harmonics,  n  =  1, 3, 5, ...,  19. 
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Fig.  5  The  normalized  intensity,  as  a  function  of  normalized  frequency,  u>/47ou;o,  and 
angle,  7o0,  of  the  radiation  scattered  by  a  relativistic  electron  (70  =  5)  from  a 
counterpropagating,  circularly  polarized  laser  pulse  (iVo  =  7,  ao  =  1-0)  for  the. 
first  three  harmonics. 


Fig.  6  The  normalized  intensity  at  the 

angle,  6,  of  the  radiation  scattered  by  a 
polarized  laser  pulse  {No  =  7,  ao~  1-0) 
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Fig.  8  The  peak  intensity  of  each  harmonic  in  the  transverse  direction  {6  =  7r/2)  versus 
normalized  frequency,  ci;/u;o,  for  a  circularly  polarized  laser  scattering  from  a 
dense  plasma  electron.  The  cases  oo  =  4  and  ao  =  6  are  shown.  The  arrows 
indicate  the  approximate  critical  harmonic  nmnber,  nc  Ug. 
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Fig.  9  The  peak  intensity  of  the  odd  harmonics  on  axis  (0  =  0)  versus  normalized  fre¬ 
quency,  a;/47oCi;o,  for  a  linearly  polarized  laser  scattering  from  a  counterstreaming 
relativistic  electron  (70  =  5).  The  cases  oq  =  4  and  no  =  6  are  shown.  The  arrows 
indicate  the  approximate  critical  harmonic  number,  Uc  3ao/4. 
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